Field Theory of Mesoscopic Fluctuations in 
Superconductor/Normal-Metal Systems 



Alexander Altland"^, B. D. Simons^ and D. Taras-Semchuk§ 
t Institut fur Theoretische Physik, Universitdt zu Koln, Ziilpicher Strasse 77, 50937 Koln, 

Germany 

= Cavendish Laboratory, Madingley Road, Cambridge CBS OHE, UK 

(February 7, 2008) 



Abstract 



Thermodynamic and transport properties of mesoscopic conductors are 
strongly influenced by the proximity of a superconductor: An interplay be- 
tween the large scale quantum coherent wave functions in the normal meso- 
scopic and the superconducting region, respectively, leads to unusual mecha- 
nisms of quantum interference. These manifest themselves in both the mean 
and the mesoscopic fluctuation behaviour of superconductor/normal-metal 
(SN) hybrid systems being strikingly different from those of conventional 
mesoscopic systems. After reviewing some established theories of SN-quantum 
interference phenomena, we introduce a new approach to the analysis of SN- 
mesoscopic physics. Essentially, our formalism represents a unification of the 
quasi-classical formalism for describing mean properties of SN-systems on the 
one hand, with more recent field theories of mesoscopic fluctuations on the 
other hand. Thus, by its very construction, the new approach is capable of 
exploring both averaged and fluctuation properties of SN-systems on the same 
microscopic footing. As an example, the method is applied to the study of 
various characteristics of the single particle spectrum of SNS-structures. 
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I. INTRODUCTION 



Physical properties of both superconductors and mesoscopic normal metals are governed 
by mechanisms of macroscopic quantum coherence. Their interplay in SN-systems, i.e. hy- 
brid systems comprised of a superconductor adjacent to a mesoscopic normal metal, gives 
rise to qualitatively new phenomena (see Ref. [p]] for a review): Aspects of the supercon- 
ducting characteristics are imparted on the behaviour of electrons in the normal region. 
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This phenomenon, known as the "proximity effect", leads to both the mean (disorder av- 
eraged) properties of SN-systems being substantially different from those of normal metals 
and various types of mesoscopic fluctuations. 

Although these two classes of phenomena are rooted in the same fundamental physical 
mechanism - a tendency towards the formation of Cooper pairs in the normal metal region 
of a SN-system - there are also major differences. Even so, it is notable that more than 
two decades passed between the first analyses of the manifestations of the proximity effect 
in the mean properties of SN-systems and the mere observation that the same effect may 
also be exhibited in mesoscopic fluctuations. The intimate connection between mean and 
fluctuation manifestations of the proximity effect will be discussed in some detail below. At 
this stage we simply itemize some basic proximity effect induced phenomena - of both mean 
and fluctuation type - and briefly comment on the theoretical approaches that have been 
applied to their analysis. 

Mean properties of SN-systems: Superconductors strongly modify the physical proper- 
ties of adjacent normal metals. For example, the proximity of a superconducting condensate 
tends to induce singular behaviour in the normal metal density of states (DoS). The prop- 
erties of these singularities depend on both the coupling to the superconductor and purely 
intrinsic characteristics of the normal mesoscopic component. This complex behaviour indi- 
cates that we are confronted with an interplay between mechanisms of quantum coherence 
in the superconductor and in the mesoscopic N-region. DoS singularities are only one of 
many more examples of manifestations of the proximity effect. For example, supercurrents 
may flow through normal metal regions and the conductance of the normal metal may vary 
with the phase of the superconducting order parameter Indeed the conductance may, 

quite counterintuitively, increase as a function of the impurity concentration, through the 
phenomenon known as reflectionless tunneling 

All these phenomena share the common drawback that they are exceedingly difficult to 
describe within conventional perturbative techniques of condensed matter physics. Broadly 
speaking, the reason for these problems is that the conventional 'reference point' of per- 
turbative approaches to dirty metals, i.e. a filled Fermi sea of electrons with an essentially 
structureless dispersion relation, represents a poor starting point to the description of SN- 
hybrids. In fact, the proximity of a superconductor leads to a strong modification of the 
states in the vicinity of the Fermi-surface, which implies that it is difficult to perturbatively 
interpolate between the conventional weakly disordered metal limit and the true state of the 
N-component of a SN-system. 

In the late sixties, Eilenberger |J also introduced a novel approach to the description of 
bulk superconductors which subsequently turned out to be extremely successful in the anal- 
ysis of SN-systems. Essentially this so-called quasi-classical approach provided a controlled 
coarse-graining procedure by which the Gorkov equation for the microscopic Green function 
of superconductors could be drastically simplified. Since on the one hand the Green func- 
tion contains all the information that is needed to describe SN-phenomena, whilst on the 
other hand - for the reasons indicated above - its computation in proximity effect influenced 
environments is in general tremendously difficult, it is clear that Eilenberger's method repre- 
sented a breakthrough. Extending Eilenberger's work, Usadel later derived a non-linear, 
diffusion-type equation for the quasi-classical Green function of dirty metals (for the precise 
definition of the term 'dirty', see below), the so-called Usadel equation. Based largely on 
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the pioneering work of Eilenberger and Usadel 0, a powerful array of quasi-classical 
methods has since been developed to study the mean properties of SN-systems. 

In view of what has been said above about the difficulties encountered in perturbative ap- 
proaches, it is instructive to re-interpret the solution of the quasi-classical equations in terms 
of the language of conventional diagrammatic perturbation theory. Referring for details to 
later sections, we here merely notice that the quasi-classical Green functions actually repre- 
sent high-order summations of quantum interference processes caused by multiple impurity 
scattering. More precisely, the solution of the Usadel equation generally sums up infinitely 
many so-called diffuson diagrams, the fundamental building blocks of the perturbative ap- 
proach to dirty metals. Whereas in normal metals high order interference contributions to 
physical observables are usually small in powers of the parameter g^ 1 [g 3> 1 is the dimen- 
sionless conductance of a weakly disordered metal), here they represent the leading order 
contribution even to a single Green function. In passing we note that, very much as the 
conductance of normal metals may be expanded in terms of the weak localization correc- 
tions, disorder averaged properties of SN-systems can be systematically expanded beyond 
the leading quasi-classical approximation in powers of g^ 1 . We will come back to this issue 
below. Having made these observations one may anticipate that the tendency to strong 
quantum interference in SN-systems not only affects their mean properties but also leads to 
the appearance of unusual mesoscopic fluctuation behaviour. 

Mesoscopic fluctuations: It has been shown both experimentally |||9|] and theoretically 
PJ^ Jl^JTl]Ji^ - [16| , p!8| , p!9| that mesoscopic fluctuations in SN-systems not only tend to be larger 
than in the pure N-case, but also can be of qualitatively different physical origin. After what 
was said above it should be no surprise that the pronounced tendency to exhibit fluctuations 
again finds its origin in an interplay between standard mechanisms of mesoscopic quantum 
coherence and the proximity effect. 

The list of examples of such SN-specific mesoscopic fluctuation phenomena includes 

• As in N-systems, sample-specific fluctuations of the conductance are universal. How- 
ever, due to proximity effect induced coherence mechanisms, the fluctuations are gen- 
erally larger than in normal systems |||§ . 

• Besides the normal current, the critical Josephson current, J c , through SNS-junctions 
exhibits fluctuations 0, which become universal in the limit of a short (L <C 



junction pOfl : 



(51?) 



j (eE c )\ L » e, 
\(eA) 2 , L«£. 



Here £ = (D/A) 1 / 2 is the coherence length of dirty superconductors, A the order 
parameter, D the diffusion constant, E c = D/L 2 the Thouless energy and L the 
system size (note that we set % = 1 throughout). 

• Such fluctuations of the supercurrent are relatively robust |J: for instance, a relatively 
strong magnetic field will exponentially suppress the average supercurrent, but reduces 
the variance of the supercurrent fluctuations by only a factor of 2. 
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The single particle spectrum in the vicinity of the DoS anomaly exhibits characteristic 
types of statistics [|l0i|T9| . 



• Novel types of universal spectral fluctuations appear JT0| . 

As compared to the mean properties of SN-systems, the physics of fluctuation phenomena 
is less well understood. Firstly, the quasiclassical approach is not tailored to an analysis 
of fluctuations. Although, to compute fluctuations, one needs to average products of Green 
functions over disorder, the quasiclassical equations are derived for single disorder averaged 
Green functions and can, to the best of our knowledge, not be extended to the computation 
of higher order cumulantsQ. Secondly, diagrammatic methods, for reasons similar to those 
outlined above, are ruled out in cases where the proximity effect is fully established. 

Important progress has been made by extending the scattering formulation of trans- 
port in N-mesoscopic systems to the SN-case [^,[11],0,^3|] . This approach made possible 



an efficient calculation of both fluctuation and weak localization contributions to various 
global transport properties of SN-systems |]J . Unlike the quasiclassical formalism, however, 
the transfer matrix approach is not microscopic. Instead, the different components of an 
SN-system are treated as black boxes which are described in terms of phenomenological 
stochastic scattering matrices. This approach, whilst extremely powerful in the analysis of 
global transport features (the conductance say), cannot address problems that necessitate 
a local and truly microscopic description. For example, it is not suitable for the calculation 
of spectral fluctuations, both global and local, the analysis of local currents, and so on. 

The purpose of this paper is to introduce a theoretical approach to the study of SN- 
systems which essentially represents a unification of the above quasiclassical concepts with 
more recent field theoretical methods developed to study N-mesoscopic fluctuations. As a 
result we will obtain a modelling of SN-systems that treats mean and fluctuation manifesta- 
tions of the proximity effect on the same footing, thereby revealing their common physical 
origin. This work represents the development of ideas that we have originally presented in 



a short letter 12 



Our starting point will be a connection, recently identified, between quasiclassical equa- 
tions for Green functions on the one hand and supersymmetric nonlinear a-models on the 
other. As was shown by Muzykhantskii and Khmelnitskii [E|J], the former can be regarded 
as the classical equations of motions of the latter. In other words, the a-model formulation 
has been shown to provide a variational principle associated to quasiclassics. So far, these 
connections have not been exploited within their natural context, superconductivity. To 
fill this gap, we will demonstrate here that by embedding concepts of quasiclassics into a 
field theoretical framework, one obtains a flexible and fairly general theoretical tool to the 
analysis of SN-systems. In particular, it will be straightforward to extend the quasiclassical 
equations so as to account for the consequences of time-reversal symmetry, the connections 



Instead of deriving equations for the Green functions themselves one may attempt to set up a 
quasiclassical approximation for their generating functional (U.Eckern, private communication). 
Whether or not such an approach has been realized and/or made working in the concrete analysis 
of fluctuation phenomena is unknown to us. 
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to perturbative diagrammatic approaches will become clear, and - most importantly - the 
effective action approach may be straightforwardly extended to the computation of meso- 
scopic fluctuations. In doing so, it will become clear in which way both the characteristic 
features of the Usadel Green function and SN-mesoscopic fluctuations originate in the same 
basic mechanisms of quantum interference. 

In this paper the emphasis will be on the construction of the approach, that is, most of its 
applications will be deferred to forthcoming publications. However, in order to demonstrate 
the practical use of the formalism we will consider at least one important representative of 
mesoscopic fluctuation phenomena, namely fluctuations in the quasi-particle spectrum, in 
some detail: The DoS of N-mesoscopic systems exhibits quantum fluctuations around its 
disorder averaged mean value which may be described in terms of various types of universal 
statistics. The analogous question for SN-systems - What types of statistics govern the 
disorder induced fluctuation behaviour of the proximity effect influenced DoS? - has not 
been answered so far. Below we will show the emergence of some kind of modified Wigner 
Dyson statistics ||25|| , within the newly constructed formalism. A concise presentation of 



both the field theory and its application to SNS-spectral statistics is contained in Ref. ||Iz 

The organization of the paper is as follows. In section [H]> we review the basic microscopic 
mechanism responsible for SN-quantum interference phenomena. In section [HI], we discuss 
the quasiclassical approach to the computation of single particle Green functions. In section 
|PV| , we briefly review the diagrammatic and statistical scattering approaches as the only 
methods so far developed to compute mesoscopic fluctuations. In the central sections |V|- 
VI l[ we introduce the aforementioned field theoretical framework. In section [V], we derive the 
effective action for a diffusive SN-structure in the form of a supersymmetric nonlinear sigma- 
model. In section [VI], we obtain the saddle-point equations of the action and examine their 
solution for some simple geometries. As mentioned above, these saddle-point equations, 
obtained by a stationary phase analysis of the effective action, are none other than the 
quasiclassical equations of motion. In section |V1J| , we address the central issue of this 
paper, the behaviour of fluctuations around the saddle-point solutions. The action displays 
a spontaneous breaking of symmetry, whose massless, or Goldstone modes are the diffusion 
modes of the system. The interaction of the diffusion modes is incorporated naturally within 
this formalism, despite their strong modification due to the proximity effect, and leads to 
mesoscopic fluctuations. We calculate in this section the renormalization of the spectrum 
of a quasi-lD SNS junction due to such fluctuations. We also demonstrate the spectral 
statistics of the SN structure to be described at low energies by a modified version of a 
universal Wigner-Dyson, or random matrix theory. The field theoretic formalism will also 
allow us to examine the onset at higher energies of non-universal corrections which serve to 
destroy the correlations described by such a universal model. In section |VIII| , we conclude 
with a discussion. 



II. ANDREEV REFLECTION AND THE PROXIMITY EFFECT 

Consider a normal metal at mesoscopic length scales, that is, scales much less than 
both Lip and Lt, where L v is the dephasing length due to electron-electron interactions 
and Lt = (T/D) 1 ^ 2 sets the scale at which the quantum mechanical coherence is cut off 
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by thermal smearing effects. The interest of such mesoscopic materials stems from the 
fact that their physical behaviour is strongly influenced by effects of large scale quantum 
interference. Such effects manifest themselves in both a variety of fluctuation phenomena 
and (non-stochastic) quantum corrections to physical observables. 

At the same time, the physics of bulk superconductors is also determined by mechanisms 
of macroscopic quantum coherence. For example, the Cooper pairs forming a superconduct- 
ing condensate represent two-electron states whose phase coherence extends over a (possibly 
macroscopic) scale set by the superconducting coherence length. 

Given that the physics of both mesoscopic metals and bulk superconductors is influenced 
by quantum coherence, it is appropriate to expect that novel interference mechanisms arise 
when two systems of this type are combined. This is indeed what happens and has led 
to the continued interest in the physics of SN-hybrid systems. A key piece of information 
required for the understanding of large scale manifestations of SN-quantum coherence is the 
manner in which normal metals and superconductors exchange quantum phase information 
on a microscopic level. The basic coupling mechanism between a superconductor and normal 
metal is a form of interface scattering, known as Andreev reflection [^,^]. In this process, 
depicted in figure ||](a), an electron at an energy below the superconducting gap, A, strikes 
the SN-interface. Due to its low excitation energy it represents a forbidden quasi-particle 
state and is unable to enter the S-region. Instead, however, it may be Andreev reflected 
off the boundary as a hole. As a result two excess charges are left at the interface which 
disappear into the superconducting condensate as a Cooper pair. 



7 



FIGURES 
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FIG. 1. (a) Andreev reflection of an electron at an SN interface, and (b) a typical pair of 
Feynman paths that lead to a non-zero value of (ip^ip^). 

The detailed physics of Andreev scattering and its consequences for SN structures has 
been reviewed extensively in the literature (see e.g. P|, |2E| , pS| -pT|] ) . Here we merely summarize 
some of its essential features that will be of importance throughout: 

• As opposed to ordinary specular reflection, Andreev-reflection represents a process of 
'retro- reflection'. More precisely, apart from a slight angular mismatch proportional 
to the excitation energy, e, of the electron above the Fermi energy, €f, the hole is 
reflected back along the trajectory of the incoming electron. 

• An electron with excitation energy e is scattered into a hole with energy — e. 

• The hole acquires a scattering phase 7r/2 — ip, where ip is the phase of the supercon- 
ducting order parameter at the interface. 

An important consequence of the existence of the Andreev scattering mechanism is the 
formation of a Cooper-pair amplitude ((ip^ (x)ip^ (x))) in the normal metal region. Here 
((...)) not only represents the quantum mechanical expectation value but also a disorder 
average. The creation of an average local pairing amplitude can be heuristically understood 
from a simple semiclassical consideration: Consider the creation of an electron somewhere at 
a point x inside a disordered metal adjacent to a superconductor (see fig. [!](&)). Due to the 
presence of disorder, the electron will propagate diffusively and may eventually strike the 
SN-interface and be Andreev reflected. In general the newly created hole may now diffuse 
along its own path. However, a particularly interesting situation arises if the hole happens 
to propagate along the path of the incoming electron back to the point of creation. As a 
result we obtain a non-vanishing pairing field amplitude {{ip 11 (x)^ (x))) . The point is that 
during their propagation through the disordered background both the incoming electron 
and the outgoing hole accumulate a quantum mechanical scattering phase which depends 
sensitively on microscopic details of the disorder. However, owing to the fact that the two 
particles propagate along the same path these phases cancel each other to a large extent. 
(For an excitation energy e = the cancellation is, in fact, perfect. For non-vanishing e one 
obtains a phase mismatch ~ L p e/vp oc eL 2 /D, where L p is the length of the scattering path, 
D the diffusion constant, Vp the Fermi velocity, L the separation of x from the interface 
and we have used the fact that for diffusive motion L p /vp oc L 2 / D.) Of course, while more 
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generic path pairs, where the electron and hole follow different paths, also contribute to 
{^y{x)^v (x)), their contributions vanish upon disorder averaging due to their strong phase 
dependence. 

The non-vanishing of (x)V> (x))) is the basic content of the proximity effect. Besides 
its resilience against disorder, the pairing field amplitude possesses a number of important 
features - all of which are related to the phase argument above - that will be of importance 
for all that follows: 

• ((i/j^ (x)^ (x))) varies weakly as a function of x. More precisely, it does not fluctuate 
on atomic scales but rather on scales set by (D/e) 1 ^ 2 . 

• ((ip^ (x)^ (x))) decays exponentially as a function of eL 2 /D. If either T or the inverse 
dephasing time r" 1 exceed e, the decay rate is set by these energy scales. 

• Quantitative expressions for the diffusive pairs of quantum paths entering the physics 
of the proximity effect are provided by so-called diffuson modes. Their meaning in the 
present context will become clear below. 

• The pairing field amplitude depends on the phases of the order parameters of the 
adjacent superconductors. If only a single superconducting terminal with constant 
phase, (p, is present, the phase dependence is simply ~ exp(i(p). In this case the phase 
is inessential and can be eliminated by means of a global gauge transformation. More 
interesting situations arise when more than one superconductor are present, in which 
case the phase sensitivity of the pairing amplitude provides the mechanism for the 
stationary Josephson effect. 

The non-vanishing of the pairing field amplitude heavily influences the properties of 
the normal metal components of SN-systems. Widely known examples of proximity effect 
induced phenomena are the DC and AC Josephson effect, which allow the possibility of 
supercurrent flows through SNS-sandwiches. Another important phenomenon is the depen- 
dence of the N-conductance on the phases of adjacent superconductors - again triggered 
by the phase sensitivity of the proximity amplitude @J5^]. However, as mentioned in the 
introduction, the emphasis in this paper will be on a study of the influence of the proximity 
effect on the single particle spectrum. 



A. Single-Particle Spectrum 

To understand the basic connection between the proximity effect and the single particle 
spectrum, let us begin by considering the simple geometry of an SNS-sandwich, shown in 
fig. 0, where the N-layer is of width L, of otherwise infinite extent and clean. This system was 
first considered by Andreev [2B| who applied scattering theory to the electron wavefunction 



to show that the spectrum in the N region, for trajectories at a fixed angle to the interface, 
is discrete below the superconducting gap. The 'Andreev levels' correspond to bound states 
with energies e given by the quantization rule, 



(eh , ,A(p\ VA 2 - e 2 

tal1 \K\ + Sgn{Vx) ~ ) = 6 ' (1) 
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where Aip the phase difference across the junction and v x the component of the electron 
velocity normal to the interface. Eq. (|l|) then indicates an average Andreev level spacing 
of the order of the inverse flight time, \v x \/L, across the normal region. A non-zero phase 
difference, Aip, leads to a shift in the levels so as to produce two separate branches of the 
spectrum, for electrons and holes respectively. 




FIG. 2. The geometry of the SNS junction. 



To determine the observed spectrum, it is necessary to sum the pole contributions to 
the DoS from bound states which arise from all possible velocity directions, according to 
Eq. (|I|). Note that non-zero contributions survive down to arbitrarily small energies due 
to trajectories travelling close to parallel to the interface. The energy dependence of the 
weights of these poles is so as to produce a total DoS which is linear in energy e, at energies 
e < A. 

The introduction of a finite concentration of impurities leads, upon disorder averaging, to 
a smearing of the formerly sharp pole structure. A less obvious outcome of this process is the 
appearance of a sharp cut-off in the spectrum, the 'minigap', below which the DoS vanishes 
entirely. The minigap, E g , is smaller than the superconducting gap, A, and depends on Aip, 
attaining its maximum for Aip = and shrinking to zero as Aip approaches it. In general 



Eg depends on Aip in a non-sinusoidal fashion [35], although in the limit of a short, diffusive 



junction, L <C £, the dependence becomes sinusoidal, E g = Acos(Ay?/2) p6| . 

In the general case, the formation of a minigap in the metallic DoS represents a highly 
non-trivial phenomenon. For example, it has been shown [37 that the presence or absence 
of a gap depends on the classical dynamical features of the metallic probe contacting the 
superconductor: For samples with integrable classical dynamics (such as the cubic system 
described above) there is no gap but rather a DoS that vanishes linearly at the Fermi energy. 
The existence of states all the way down to zero energy may be understood by means of 
Bohr-Sommerfeld quantization arguments. By contrast, in systems with chaotic classical 
dynamics, a gap opens whose magnitude may depend on both the coupling strength to the 
superconductor and the intrinsic classical transport time through the metal region. 



B. Failure of Semiclassics: Quantum Diffraction 



The prediction of a minigap provides a useful test for any analytic approach to the 
physics of SN-systems: for example, an application of standard approximation schemes to 
semiclassical formulae for the the DoS fails to predict correctly a gap. 
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To explain this point, let us consider the simple case of a diffusive metallic cube of linear 
extension L attached to a superconductor. In this case the DoS gap is of width ~ E c . 
In order to embed the previous heuristic path arguments regarding the proximity effect 
into a quantitative calculation of the proximity effect influenced DoS one might apply the 



semiclassical Gutzwiller trace formula ||38|| , a powerful computational tool often used in the 
analysis of the DoS and correlations thereof. The Gutzwiller trace formula essentially states 
that the averaged DoS can be obtained as a sum over all periodic orbits with vanishing, or 



at least a disorder insensitive action 3£ . Small action contributions to trace formulae are 



usually computed from the so-called diagonal approximation. In the present context the 
diagonal approximation would amount to counting Feynman paths such as the one depicted 
in fig. g|(a). (The smallness of the action of these pairs of paths follows from the fact that if 
the action of the electronic segment of the path (solid line) has an energy dependent action 
S(e), the action of the hole segment (dashed line) will be —S(—e). The two contributions 
nearly cancel each other.) 

The problem with the standard diagonal approximation is that it fails to predict a gap 
in the DoS, even if paths with multiple Andreev scattering are taken into account. In order 
to understand this failure we first have to notice that small action path-pairs exist which do 
not fall into the scope of the diagonal approximation (by which we mean that they cannot 
be obtained as a superposition of two identical segments, one electron- one hole-like). A 
common feature of these 'non-diagonal' path configurations is that they contain 'junction 
points' were the paths of electrons and holes split (cf. fig. 0(6)). In order to understand the 
existence of these splittings one has to keep in mind that the paths entering the semiclassical 
picture do not correspond to rigorously defined solutions of classical equations of motion but 
should rather be thought of as objects that are smeared out (in configuration space) over 
scales comparable with the Fermi wavelength jHJ. As a result two classically 'identical' 
paths may split and recombine at some later stage, a process which is not accounted for 
by the diagonal approximation. This splitting, as it is caused by the wave nature of the 
electrons, is sometimes referred to as a quantum diffraction phenomenon. 

Note that the junctions appearing in fig. |3](6) are reminiscent of similar processes needed 
to generate weak localization corrections to the conductance of normal metals jllj . However, 
whereas weak localization corrections represent a correction of 0{g~ x ) to the classical con- 
ductance (g ^> 1 is the dimensionless conductance), the diffraction corrections appearing in 
the present context can by no means be regarded as small. In fact they are as important as 
the leading order diagonal contributions which implies that processes with up to an infinite 
number of 'junction points' have to be taken into account. This fact not only explains the 
failure of the diagonal approximation but also the difficulties encountered in diagrammatic 
analyses of the proximity effect. The point is that each of the 'legs' appearing in fig. $(b) 
represents a Cooperon. The perturbative summation of infinitely many Cooperons repre- 
sents a difficult problem, in particular in cases where the sytem is truly extended in the 
sense that it cannot be treated within an ergodic or zero mode approximation (for a per- 
turbative analysis of the zero- mode scenario, see Brouwer et al., [5T]]). Fortunately there is 



an alternative approach, the quasiclassical method reviewed below, which provides a highly 
efficient tool for the effective summation of all interference corrections contributing to the 
DoS and other physical observables. 
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(a) (b) 
FIG. 3. Trajectories that are (a) included and (b) not included within a semiclassical treatment. 

III. QUASICLASSICS 

In order to prepare the discussion of quasiclassics we first need to introduce its micro- 
scopic basis in general and the Gorkov Green function in particular. 

A. The Gorkov Equations 

As long as interaction effects are neglected^ the complete information on any SN-system 
is encoded in its single particle Gorkov Green function, Q r,& . The latter is defined by the 
Gorkov equations |EJ, whose matrix representation reads 

/i + e ± -^(p-fA( ri )) 2 -V( ri ) A(n) 

-A(rr)* /i - e± - ^ (p + f A( ri )) 2 - Vfa) 

= M d ( ri -r 2 ), (2) 

where 

Ft r,a G t r ' a J' ( 3 ) 

Here G r,a and F r,a represent the normal and anomalous Green function, respectively, A(r) 
is the superconducting order parameter (which, in principle, has to be determined self- 
consistently), e± = e ± iO, and V is the impurity potential. 

The equation above may be represented in a more convenient way by introducing Pauli- 
matrices, <jf h , operating in the two-component particle/hole spacef]. Separating the order 



2 In this paper the - important - role Coulomb interactions may play in SN-physics will not be 
discussed. 

3 Recall the general definition of the Pauli matrices, 



e? r ' a (ri,r 2 ) 
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parameter into its modulus, |A(r)|, and phase, y?(r), Eq. (|2|) takes the form 

g^(r 1 ,r 2 )=5 d (v 1 -r 2 ), (4) 

where A = af |A|e^ ri ) ff 3 . The presence of an impurity potential, V(r), makes the solu- 
tion of Eq. (|4]) difficult. However, as pointed out by Eilenberger |J and Larkin and Ovchin- 
nikov , a crucial simplification applies in the case where the wavelength of the electrons, 
Xp, is small as compared to the characteristic scales over which the order parameter A(r) 
and vector potential A(r) vary. Under this condition one can resort to the 'quasiclassical 
approximation' . 



e F -±-(p- ^A(rK h ) 2 - V(t) + (A(n) + 6 ± ) al 



B. Quasiclassical Approximation 

The starting point of the quasiclassical approach is the observation that the spatial 
structure of the Green function is comprised of rapid oscillations, over a spatial scale of 
the Fermi wavelength, modulated by a slowly fluctuating background over longer scales. A 
quasiclassical analysis of the proximity effect involves an averaging over the rapid variations 
of the Green function. At the same time, sufficient information is retained within the slower 
modes of the Green function to provide a useful approximation to the full consequences of 
the proximity effect. The advantage is a great simplification of the corresponding kinetic 
equations. As the derivation of the quasiclassical equations has been reviewed extensively 
in the literature (e.g. [H,[3T,49]]), we restrict ourselves here to a brief summary of the main 



results of the approach. 

The quasiclassical, or Eilenberger, Green function, g r ' a (n, r), is obtained from the Gorkov 
Green function by a) a Wigner transform, b) an impurity average and c) an integral over 
the kinetic energy variable. The precise definition reads 

<7 r ' a (n, v)= l -jdij d( ri - r 2 )£ r < a (r 1; r 2 ) expHp • (n - r 2 )), (5) 

where r = (rx + r^)^, £ p = Vp(p—pF), n = p/p and ^Op = mvp is the Fermi momentum. The 
application of this approximation to the Gorkov equation, Eq. (f|) leads to the 'Eilenberger 
equation' B: 



v F • V r # r ' a (n, r) 



af(e ± + A(r)) + - <^ a (n', r)) n , , r) , (6) 



where r is the elastic scattering time due to impurities. This equation essentially represents 
an expansion to leading order in the ratio of to the scale of spatial variation of the 
slow modes of the Gorkov Green function. It can be shown that the Eilenberger Green 
function obeys the nonlinear normalization condition (see, for example, the discussions in 
Refs. |7l| ) 

<f> a (n,r) 2 = ll. (7) 

Eq.(^|) represents an equation of Boltzmann-type which is much simpler than the original 
Gorkov equation, but may still be difficult to solve in general. However, significant further 
simplifications are possible in the 'dirty limit'. 
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C. Dirty Limit 



The 'dirty limit' is specified by the conditions i = vpT -C £ (implying that the dominant 
transport mechanism is diffusion) and e < r _1 (implying that 'time scales' e _1 much longer 
than the scattering time are explored). Under these conditions, the dependence of the Green 
function on the angular direction (represented by n) is weak and one may expand in its first 
two spherical harmonics: 

^ a (v F) r)=^ a (r) + n-g r 1 ' a (r) + ..., (8) 

where go' a ( r ) ^ n ' gi' a ( r )- A systematic expansion of the Eilenberger equation in terms of 
gi then leads to a nonlinear and second-order equation for the isotropic component, 

£V(^' a (r) V«£ a (r)) + ^ 3 ph (e ± + A(r)), <£ a «] = 0, g^ivf = 1, (9) 

known as the 'Usadel equation' 0. In order to specify a solution, one has to supplement the 
equation with appropriate boundary conditions. The analysis of the boundary behaviour of 
the equation becomes somewhat technical. For this reason a more detailed discussion of the 
boundary conditions 

v(-)god r go(-) = v(+)god r go(+), (10a) 

a(±)g d r9o (±) = ^[go(+),9o(-)], T « 1, (10b) 
0o(+) = </oH, (10c) 

has been made the subject of appendix [A|. Here T e [0, 1] is a measure for the transparency 
of the S/N interface, <7o(i) denotes the Green functions infinitesimally to the left respectively 
right of the junction, and Gt (cf. Eq. ( |A5| ) ) is the tunnel conductance of the interface. 



D. Solution of the Usadel Equation 

Solutions of the Usadel equation with appropriate boundary conditions have been derived 
for a vast number of geometries. At the same time a systematic and general solution scheme, 
based on an effective circuit theory, has been constructed by Nazarov [52]. Furthermore, 



a number of quasiclassical predictions seem to be borne out well experimentally (see e.g. 
|34] , |53| -p7| . In the field theoretic context introduced below, the Usadel equation and its 
boundary condition will reappear on the level of the mean field analysis in section [VT]. In 
addition explicit solutions for some simple geometries will be discussed in that section. 

It is instructive to consider a diagrammatic reinterpretation of the Usadel solution in 
terms of a summation over real-space trajectories. Such a decomposition may be achieved 
by taking the solution of Eq. (^) in a normal region, with the SN interface (at x — 0, say) 
represented (for energies e < A) by a boundary, or source, term AS(x). An expansion of 
the solution in powers of A corresponds to a series of diffusive trajectories which include 
successive numbers of Andreev reflections at the interface. Fig. [5] (b) provides an illustration 
of a trajectory with four reflections, whilst 'starfish' trajectories with arbitrary numbers 
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of reflections are clearly possible. The reproduction of the Usadel solution requires the 
summation of the full set of trajectories, corresponding to diagrams to all orders. We see 
that the inherent difficulties of a diagrammatic treatment of the proximity effect, as discussed 
earlier, extend even to the relatively simple task of reproducing quasiclassics. 



IV. BEYOND QUASICLASSICS 

The quasiclassical approach allows for the efficient calculation of a wide spectrum of 
physical observables. In general, any observable that may be expressed in terms of a single 
disorder averaged Green function is a candidate for quasiclassical analysis. Note that, by 
extending the formalism so as to include Keldysh-Green functions [[58 1, observables that 



are commonly expressed in terms of two-particle Green functions also become accessible. 
Important examples are the conductance and other transport quantities. However, there 
are important classes of observables which do not fall into the above category, thereby 
falling beyond the scope of quasiclassics. The list of inaccessible quantities may be grouped 
roughly into four different categories: 

• Physical observables which, by definition, are not expressible in terms of single particle 
Green functions. An example is given by the magnetic field dependence of the London 
penetration depth for a bulk superconductor, as studied by Larkin and Ovchinnikov 
|43|j . Its analysis requires the computation of the average of four momentum opera- 
tors, (p(0)p(t)p(0)p(t)) , a quantity that involves two- rather than one-particle Green 
functions. 

• Higher order quantum interference corrections to single particle Green-functions. Very 
much like weak localization corrections of 0(g~ l ) to the classical Drude conductance, 
the quasiclassical Green function represents the leading order term of a series expansion 
in powers of g~ x . The next to leading order contributions become important in cases 
where one is interested in quantum corrections of weak localization type or strong 
localization effects. 

• The quasiclassical approach (in its extension to include a Keldysh component) does 
not account for the corrections to two-particle Green functions due to the interference 
of mutually time-reversed trajectories. An example property that is affected in this 
way is the conductivity, as we see below. 

• Most importantly, the quasiclassical approach does not allow for the study of meso- 
scopic fluctuation phenomena. The analysis of fluctuations requires the computation 
of disorder averages of two or more Green functions. Due to the impurity induced in- 
terference between different Green functions, quasiclassical techniques are inapplicable 
to these problems. 

Given that these classes of problems cannot be addressed within quasiclassics, it becomes 
necessary to seek some alternative approach. Here, we briefly review two of perhaps the most 
important theoretical tools currently established, namely, diagrammatics and the scattering 
matrix approach. 
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A. Perturbative Diagrammatic Methods 



Microscopic diagrammatic methods have been applied to the study of various SN- 
phenomena. The list of diagrammatic analyses includes computations of Josephson cur- 
rent fluctuations through SNS junctions ||, investigations of the phase sensitivity of the 
N-conductance ||, computations of universal conductance fluctuations of SN-systems 
and more. 
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FIG. 4. (a) shows the trajectories whose interference leads to the first weak localization cor- 
rection to the conductivity. (6) shows the corresponding diagrams, where the curvy line is the 
cooperon and the triangles represent Andreev reflection processes. 

The basic building block of diagrammatic analyses are the diffusion modes both of dif- 
fuson and cooperon type. What makes these modes different from their counterparts in 
pure N-systems is that they now include Andreev scattering events, as represented in fig. |j. 
In the presence of the proximity effect, any of the Andreev scattering vertices appearing 
in these modes is in turn to be renormalized by further diffusion modes, as indicated in 
fig. |^. Note that these diagrams are the formal representation of what in real space are the 
'legs' of the starfish-like structures appearing in fig. || The problem with the diagrammatic 
approach is that in situations where the proximity effect is fully established, the Andreev 
vertices renormalize heavily, i.e. one has to sum self-consistently nested series of the dia- 
grams depicted in fig. [5[ Another way of putting this is to say that one has to perturbatively 
reconstruct the solution of the Usadel equation, a difficult if not impossible task. In fact, 
an incomplete account of the proximity-induced renormalization processes may lead to in- 
correct results: for example, diagrammatic analyses of universal conductance fluctuations 
by Takane et al. [|ll|] failed to reproduce correctly their surprising insensitivity to external 
magnetic fields, as later demonstrated by Brouwer et al. |jl], |16|l . However, in cases where the 
proximity effect is either suppressed or of secondary importance, diagrammatic tools can be 
applied successfully to the study of SN-systems. 

In summary, it can be said that diagrammatics is applicable to the perturbative analysis 
of SN-phenomena in cases with a weakly pronounced proximity effect. As is usual with 
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diagrammatic methods, non-perturbative problems, such as localization, fine structure level 
statistics, and so on, cannot be addressed. 




FIG. 5. First order renormalization of the Andreev scattering process necessary to include the 
proximity effect. 



B. Multiple Scattering Formalism 

Scattering theory provides a powerful theoretical tool for the analysis of quantum trans- 
port through mesoscopic systems in general, and SN-systems in particular. The reason for 
the efficiency of the scattering theoretical formalism in the study of SN-systems is not only 
its relative simplicity, but also the fact that the proximity effect does not seem to lead to 
essential complications. Due to this latter advantage, scattering theory has for a long time 
existed as the only tool for computing mesoscopic fluctuations under the influence of the 
proximity effect. 

The starting points of the scattering approach are generalizations of the standard multi- 
channel Landauer formulae for N-mesoscopic systems to the SN-case. For example, in the 
case of a single N-sample attached to a superconductor the conductance may be expressed 



as [P,y.[22 



2e 2 4e 2 

GnS = ~^ _tr ( 1 ~~ Beetle + She/he) = -J^^^he^he^ 

where S ee (S e h) are the matrices describing the scattering of electrons incoming from the nor- 
mal metal to electrons (holes). In a next step, the scattering matrices are expressed in terms 
of a) the transmission matrices of the normal metal compound (which are known in terms 
of their transmission eigenvalue distribution functions P|,p9[] , and b) matrices describing the 
scattering off the superconductor. 

The scattering theoretical approach is particularly powerful if the observables of interest 
take the form of 'linear statistics', i.e. quantities, X, that can be represented as 

X = £/(T n ), 

n 

where / is some function and T n is related to the n-th normal transmission matrix eigenvalue. 
This is often but not always the case: for example, in time reversal invariant cases, the 
conductance of the above SN-system may be formulated as (see Beenakker, 



4e 2 N T 2 

^"^(rf^- (12) 

However, if time reversal invariance is broken, such a simple representation is no longer 
possible and expressions involving not only eigenvalues but also the diagonalizing matrices 
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appear. It has been shown by Brouwer et al. that, even under these more compli- 
cated circumstances, scattering theory remains applicable. What becomes necessary is to 
supplement the conventional transfer matrix approach (by which we mean the derivation 
and solution of a Fokker-Planck equation for the eigenvalues) by diagrammatic methods 
accounting for the presence of the diagonalizing matrices. Due to the complex structure 
of the diagrammatic series appearing in the SN-problem, only quantum dots (i.e. ergodic 
systems that can be modeled in terms of a scattering matrix distributed through a single 
Haar measure) rather than arrays thereof could be analyzed in this way: we again encounter 
the notorious difficulties accompanying perturbative approaches to SN-systems. 

In passing we note that in some cases SN-quantum dots, chaotic or disordered, can be 
modeled in terms of simple random matrix theory. More precisely, random matrix techniques 
become applicable if the proximity effect is suppressed. This happens if the system is subject 
to a magnetic field (of the order of a few flux quanta through the system) , or if the phases of 
the adjacent superconductors average to zeroQ In spite of the fact that the proximity effect 
is suppressed, the mechanism of Andreev scattering remains active and manifests itself in 
the SN-random matrix ensembles having symmetry properties that differ substantially from 
the standard Wigner Dyson ensembles. 

To summarize, statistical scattering theory represents a powerful tool for the analysis of 
both mean and fluctuation characteristics of global transport quantities. Clearly, observables 
belonging to this category are of outstanding importance from the experimental point of 
view. Nevertheless, problems remain for which one is interested in observables that are 
local and/or microscopically defined. This complementary class of quantities is inaccessible 
through phenomenological scattering analyses. Thus at least one alternative theoretical tool 
for the analysis of SN systems is called for. 



V. FIELD THEORY FOR SN SYSTEMS 

In the following central part of the paper we are going to introduce a novel approach to the 
analysis of SN systems which is based on field theoretical and, by construction, microscopic 
concepts. The formalism will be applicable to observables that can be expressed in terms 
of one or products of single particle Green functions. If this criterion it met, both mean 
values (including quantum corrections to quasiclassical results) and mesoscopic fluctuations 
can be computed. In a few exceptional cases, distribution functions can be obtained. It will 
quickly turn out that the formalism is intimately related to each of the approaches reviewed 
above: On the mean field level it reproduces quasiclassics, perturbative fluctuations around 
the mean field can be interpreted diagrammatically, and the connection to scattering theory 



4 The latter mechanism is rather subtle and it is not clear whether it can be realized in practice. 
The reason is that even minute phase fluctuations of 0(g^ 1 ) invalidate the applicability of random 
matrix theory. However, systems with a phase-suppressed proximity effect are realized in nature 
as vortices in type II superconductors, see e.g Caroli et al. [ SO | . Here the vortex center has a 
non- vanishing metallic DoS, an indication of proximity effect suppression. 
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is established through general parallels between cr-models and the transfer matrix approach 

n. 

Prior to embarking on any kind of detailed discussion, let us briefly outline the main 
conceptual steps of the construction of the field theory: 

1. Starting from the microscopic Gorkov Hamiltonian of an SN-system we will construct 
a generating functional for the disorder average of the product of a retarded and 
an advanced Green function. (Generalizations to products of more than two Green 
functions are straightforward.) The functional will be of a nonlinear cr-model type. 
Essentially it represents a supersymmetric extension of earlier (replicated) field the- 



oretical approaches to bulk superconductors [66]. In the present formulation, the 
order-parameter field will be imposed and not computed self-consistently (a common 
and mostly inessential simplification in the field of SN-systemsQ) . 

2. It will then turn out that the spatial inhomogeneity of the order parameter field poses a 
substantial problem: A straightforward perturbative evaluation of the field theory, by 
which we mean a perturbative expansion around any spatially homogeneous reference 
field configuration, is impossible. It goes without saying that this is nothing but the 
manifestation, in a field theoretical context, of the general perturbative difficulties 
characteristic to SN-systems. 

3. The way out will be to subject the field theory, prior to any perturbative manipula- 
tions, to a stationary phase analysis. Given what has been said under the previous 
item, it is clear that the saddle point configurations of the theory must be spatially 
inhomogeneous. More specifically it will turn out that the stationary phase equation of 
the theory is simply the Usadel equation. In other words, the quasiclassical approach 
to SN-systems turns out to be equivalent to the mean field level of the field theoretical 
formalism. 

4. We then turn to the issue of fluctuations. Broadly speaking, two qualitatively different 
types of fluctuations will be encountered: Massive fluctuations around the mean field 
(giving rise to quantum corrections to quasiclassics) and a Goldstone mode. The latter 
will induce correlations between retarded and advanced Green functions and thereby 
mesoscopic fluctuations. 

To keep the discussion of the above hierarchy of construction steps from being too abstract, 
the computation of correlations in the single particle spectrum will serve as a concrete 
example of an application of the theory. 



5 Self-consistent calculations of A(r) in mesoscopic SN-junctions, as achieved analytically by Zaikin 
pl| in the clean case and numerically by several authors (e.g. Refs. [69-f7l|]) in the dirty case show 
that A is suppressed in the S region near the interface. It may also become non-vanishing in the 
N region if an electron-electron interaction is included there. We will neglect these effects as we 
do not anticipate that they would have a significant impact on our results. 
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Before turning to the actual construction of the field theory, it should be noted that 
essential components of the machinery we are going to discuss are not original, but have been 
introduced earlier: The general supersymmetric field theoretical approach to N-mesoscopic 

H 



systems has been constructed by Efetov 



Oppermann ||66|| , and later Kravtsov and 



Oppermann ||67|| , introduced a fermion-replicated a- model description of bulk disordered 
superconductors. As far as technical aspects are concerned, the formalism we are deriving 
represents a supersymmetric extension of Oppermann's model, tailored to the description of 
spatially inhomogeneous structures. In its early stages, the construction of the model follows 
a by now absolutely standard strategy. Essentially, this requires an adaptation of Efetov's 
model to allow for the particular structure of Gorkov Green functions. For this reason our 
presentation of the early construction steps will be concise, but nevertheless self contained. 



A. Field Integral and the Ensemble Average 

As is usual in the construction of field theories of mesoscopic systems, the first construc- 
tion step is to represent products of matrix elements of Gorkov Green functions in terms of 
supersymmetric Gaussian field integrals. To keep the discussion comparatively simple, we 
focus on the case of two-point correlation functions {G T (e + u + /2)Q a (e — u + /2)) and choose, 
as a specific example, the quantity 

(tr pV (G T (e + u+/2)) tr ph , r (£ a (e - u+/2))) , (13) 

which appears in the computation of the fluctuations of the DoS, Q 

R(e, co + ) = 7-^73 (^(e + u+/2)6u(e - u+/2)). (14) 
(z/(e)) 2 

Here 8u = v — (u), u(e) = — ^Tmtr phr (<? r (e)) and tr phr denotes a trace with respect to 
both position and particle hole index. 

In order to represent objects of this kind in terms of Gaussian field integrals, we first 
introduce a 16-component vector field, ip — f0A,s,a,f( r )}j A, s, a, t = 1,2, with complex 
commuting (anticommuting) components a = 1 (a = 2). The significance of the two- valued 
indices A, s, a, t is summarized in the following table: 



index 


significance 


abbreviation 


A 


advanced/retarded 


ar 


s 


particle/hole 


ph 


a 


boson/fermion 


bf 


t 


time reversal 


tr 



6 Note that the correlation function i?(e, u+) differs in two respects from the analogous quantity 
in ./V-mesoscopic systems: (i) The mean DoS, (v(e)}, will not, in general, be constant. Hence, 
R(e,u+) may explicitly depend on the center coordinate e, a fact that can be remedied by an 
unfolding procedure (see e.g. |38|). (ii) as opposed to iV-systems, correlation functions such as 
(G r G r ) are non-trivial in the sense that they do not equal the product of averages. Both aspects 
(i) and (ii) will be commented on later in more detail. 
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Apart from the index s accounting for the 2x2 matrix structure of the Gorkov Green 
function, all other indices are standard in supersymmetric approaches to disordered systems. 
For a discussion of their significance we refer to Efetov's book We next introduce the 
action 



S 

where 



■0 



fi - — p - -Act 

2m 



P h 



CT 



tr 



V + a 



P h 

3 



A + e + — ct 3 



1>, 



A(r) = Act] 311 <g> of exp(-^(r)orf n <g> ct 



ph 



f), 



(15) 



(16) 



and the Pauli matrices crf r , crf h , CTj bf , and ct 



bf 



= 1,2,3 operate in the two dimensional 
spaces of A, a, s, t-indices respectively. The fields ip and tjj are related to one another via 

where the matrices Efj, x = ar, ph, bf,tr, are defined as [EfA = 5ut5jjr, and the indices 

i',f refer to the space x. 

As in analogous theories of iV-systems, the action, S, can be employed to represent 



correlation functions in terms of Gaussian field integrals. Specifically, the expression (|13"D 
takes the form 



(tr ph , r (g r (e + UJ+/2)) tr ph , r {Q\t - w+/2))) 
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^af®af^ / ^ 2 CT 3 bf ®CT 3 P V 2 , 



(17) 



where the indices on the ^-fields refer to the ar-space, and all other indices are summed 
over. 

After this preparation - which essentially has comprised of an extension of existing 
supersymmetric framework to account for the additional ph-structure - we may proceed in 
strict analogy to standard procedures: 

• Firstly, averaging over Gaussian-distributed disorder, with the correlation function 

1 



(V(r)V(r')) 



2nv n T 



5(r-r>), 



where v n denotes the DoS of a bulk N system, generates the quartic contribution to 
the action, 



1 



int. 



47TIVT 



Next, Sint is decoupled by introducing a 16 x 16 Hubbard-Stratonovich matrix field, 
exp(-S , int ) = J VQ exp J (^pQip - ^strQ 2 
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where 'str' denotes the supersymmetric extension of a matrix trace[|. 
• In a third step we integrate over the ■f/'-nelds to arrive at the Q-represented action, 

S[Q] = istr.hKT 1 " ^ / strQ 2 , (18) 

where 

G- 1 = A* - ^ (p - - c ^f ® ^ + ^ ® (A(r) + e + + (19) 

and str r denotes a trace extending over both internal and spatial degrees of freedom. The 
next step in the standard construction of the a- model is to subject the action S[Q] to a 
saddle point analysis. As we will see shortly, the presence of a superconductor with spatially 
inhomogeneous order parameter will necessitate substantial modifications to the standard 
mean field scenario. In order to gain some insight into the structure of the mean field 
equations, the next section will be devoted to the study of the comparatively simple case of 
a bulk disordered superconductor. However, prior to specializing the discussion, let us make 
some general remarks as to the structure of the stationary phase equations. 



Varying the action fll8|) with respect to Q generates the stationary phase equation, 



i 



Q(r) = g(r,r). (20) 

7TI/„ 



Before embarking on the explicit computation of solutions to Eq. (|20f) - which represents a 
16-dimensional matrix equation - it is convenient to elucidate further its structure. In fact, 
a striking simplification arises from the fact that the Q-independent part of the kernel Q~ x 
is diagonal in all indices save the ph-indices. 

In the standard case, that is, no superconductivity and hence no ph-indices, the diag- 
onality is complete and one may start out from an ansatz for Q which is fully diagonal. 
Exploiting the fact that the energy- difference between the Green functions, u+, is typically 
small in comparison with all other energy scales of the system, one might be tempted to 
argue that Q~ l is not only diagonal but even approximately proportional to the unit matrix 
in the internal indices. As a consequence one might assume that Q is proportional to the 
unit matrix as well. This, however, is wrong. The infinitesimal imaginary increment con- 
tained in u + gives rise to a phenomenon of spontaneous symmetry breaking in the sense that 
the saddle point solution in the retarded sector differs from the one in the advanced sector 
(see e.g. More precisely, the saddle point solution of the field theory for iV-systems 



reads^ 



7 We use the convention that str A = tr^4bb — trAg. 

8 Note that <7f r is commonly denoted by A in the literature. 
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What kind of modifications arise in the presence of superconductivity? Firstly, the kernel 
Q^ 1 is no longer fully diagonal. It contains a non-trivial matrix structure in ph-space[]. Thus, 
the simplest ansatz for a saddle point solution is diagonal in all indices save the ph-indices. 
Secondly, we may expect that, as in the N-case, the structure of the solution depends on 
the infinitesimal increments added to uj + . An inspection of Eq. fll5| ) tells us that the role 
played by the matrix af in the standard case will now be taken over by of 1 <S> 0"f r . Finally, 
the solution in the N-case, Qn = ofS was fully universal in the sense that it did not depend 
on any energy scale. Since the 'perturbations' arising in the action due to the presence of 
the superconductor - most notably the order parameter - are weak in comparison with the 
dominant energy scale, fi, it is sensible to assume that the eigenvalues of the saddle point 
solution will still be ±1. 

Starting from the comparatively simple example of a bulk superconductor we will 
next confirm these suppositions by explicit calculation. 



B. Example: A Bulk Superconductor 

In this section we assume the order parameter to be spatially constant. The resulting 
saddle point equation has previously been discussed in pBL Specifically, we assume A(r) 



|A|a2 h , corresponding to a constant gauge ip = tt/2, and assume that the vector potential 
vanishes, A = 0. Homogeneity of the gap function implies homogeneity of the saddle-point 
solution Q(r) = Q. We next introduce the ansatz 




Q = q • a** 

where <r ph = (o^ , , o^) T is a vector of Pauli matrices operating in ph-space. The com- 
ponents of the vector q are diagonal matrices which are trivial in all but the ar-space, 

(21) 

and normalized to unity, q± • q± = 11. To proceed it is convenient to adopt an elegant 
parametrisation for the Green function suggested by Eilenberger ||. First note that in 
momentum representation Q~ l can be written as 

Q-\v) = -C( P ) + iw-^ h , 

where 



9 One might argue that the off-diagonality may be removed via a global unitary transformation, 
at least in the case of a spatially homogeneous order parameter. This, however, would contravene 
the conditions enforced by analyticity on the structure of the imaginary increments contained in 
the action, which, as we saw in the N-case, play a crucial role in determining the structure of the 
solution. 
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2m ' 
q 

w = - ir, r = (e + i0a^ v )e 3 + iAe t , 

It 

and u+ has been set to zero. Using the fact that (w • <r ph )(w • cr ph ) = w • w = |w| 2 , it is a 
straightforward matter to show that 

1 l ph + sw ■ qP h 

2 ,t±i -C(P) +*s|w| 

where w = w/|w|. Performing the trace over momenta, and making use of the relations 

^ / (l^" /i )^ Vol ' Zy "/ rfC/(C) ' (23) 
where 'Vol' is the system volume, as well as 

s ' n _ J n = 0, 

^ _ C + ia | w | - \ -2m + 0(max (A, e)//i) n = 1, (24) 

the saddle point equation fl20| ) takes the simple form 



-£( r > r ) = r^r£S(p) 



7TZ/„ 7TZ/ n Vol 

•£4> q = w. 



w • er 



q = " (25a) 



which is solved by 



R=\A\ 2 -e 2 . (25b) 

To achieve the correct analyticity of q, we take the branch cut in the square root in R to be 
along the negative real axis. Introducing the parameterization 

q = (sin# s , 0, cos# s ) , (26) 

we find 



cos6L 



H > |A|, 



R\y* l-zsm(e). lei < |A| 



s ™ s = Tkm. AL,^ a ( 27 ) 



_i?|i/2 l-isgn(e), |e| < |A|. 

In the particular limit of zero order parameter, the solution collapses to a ph-diagonal one, 
q = e^of ~* Q = erf 1 <g> a| r . This result is consistent with the conventional saddle-point 
equation of the normal conductor (cf. the remarks made towards the end of the preceding 
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section): For vanishing order parameter, the particle/hole extension simply generates two 
copies of the normal Hamiltonian. The effect of a non-vanishing value of A is to induce a 
rotation of the saddle-point in the ph sector. In the extreme limit of e — > with a finite 
gap, the saddle point rotates as far as q = sgn(e)ei, or 9 S = sgn(e)7r/2. 
The local DoS can be computed from q asQ 

(i/(r)> = -^Imtr ph (^(e; r, r)> = ^ Re (str {Qaf ® <x 3 ph ® E™)) Q Q -=$ u n Re[q + (r)] 3 . (28) 

Here (■ ■ -)q denotes the functional expectation value of the field Q and the final expression 
is obtained by evaluating the functional integral at its saddle point value. For the bulk case, 
this gives a superconducting DoS, u s , of u s = u n Recos9 s , which leads by Eq. (|27D to the 
familiar BCS form, 

u 8 = e(|e| - |A|) 



(|A|2_ e 2)l/2- 

Note that below the gap, |e| < |A|, the loss of antisymmetric ar structure of q in Eq. ( |27|) 
is reflected in a zero DoS. 

We may next ask whether the saddle point solution, q, is unique. Anticipating that 
all saddle point configurations must share the eigeni>a/iie-structure of Q, a general ansatz 
probing the existence of alternative solutions reads 

Q -> TQT-\ (29) 

where T is some rotation matrix. In order to understand the structure of the resulting saddle 
point manifold, it is essential to appreciate that there are three parametrically different 
energy scales in the problem. 

• The asymptotically largest scale in the problem is the chemical potential. The existence 
of the large parameter fi/E, where E may be any other scale involved, stabilizes the 
eigenvalue structure of the matrix Q. (This follows ultimately from the structure of 



the pole integral fl24|) - see also the corresponding discussion in ||.) Thus, as long as 
we are not interested in corrections of (D(E/ fj,), configurations of the type of Eq. fl2"9"|) 
exhaust the field integration domain of interest. 

The next largest scales are A and/or e. Amongst the configurations parameterized by 
Eq. (|2~9p, there are some that are massive in these parameters and some that are not. 



10 To derive Eq. ( |28| ) one starts out from the functional representation of the local DoS, 

tr ph (S r (e;r,r)> = ±J Vty, ^e'^^af ® of ^(r). 

After the Hubbard-Stratonovich transformation the preexponential terms take the form of a func- 
tional expectation value ~ (str (Qa\ { <8> erf 1 (g> E^))q which, upon evaluation in the saddle point 
approximation leads to Eq. (EH). 
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• The smallest scale is u+. Its physical significance is that of an inverse of the time scales 
at which we are probing correlations. With regard to correlation functions in u + , field 
fluctuations that are massive in the intermediate parameters e and A are clearly of 
little if any relevance. 

After these preparatory remarks, it should be clear that we will be concerned mainly with 
fluctuation matrices T that still lead to solutions of the saddle point equation up to correc- 
tions oc uj + . An inspection of the action (|18|) tells us that such T have to fulfill the condition 
[T, <j ph ] = 0. On the other hand the matrices T must not commute with Q, as otherwise 
they would be ineffective. Taking these two facts together, we see that the most relevant 
fluctuations, T, around Q, are those that generate rotations in ar-space: For u + = any 
configurations fulfilling the above conditions again represent solutions of the saddle point 
equation. In other words, the T's generating these configurations are Goldstone modes. 

Before extending the discussion to spatially inhomogeneous problems and the impact 
of the existence of Goldstone modes, let us comment on a mathematical aspect of the 
construction of the theory. We have seen that, for A 7^ 0, the saddle point solution Q 
differed substantially from the standard saddle point erf 11 (g> crf r of the bulk metallic phase. 
This raises the question of whether the superconducting saddle point, q ■ cr ph , and erf 11 ® 
are both contained in the field manifold of the nonlinear cx-model. Clearly this question will 
be of concern as soon as we deal with SN-/i?/6noJ-systems, and, in fact, the answer is negative. 
However, it turns out that the problem can be surmounted by analytic continuation of the 
parameter space spanning the Q-field manifold. Since the discussion of the manipulations 
needed to access both saddle points is inevitably somewhat technical, it has been deferred 
to an appendix and may be skipped by readers who are not interested in details of the 
formalism. 

Our final objective will be to describe SN-systems rather than bulk superconductors. 
What makes the analysis of SN-systems technically more involved is that the action is man- 
ifestly inhomogeneous: In our comparatively coarse modelling, where the order parameter 
is externally imposed, the N-component of an SN-system will be characterized by a sud- 
den vanishing of the order parameter. Within a more accurate description, based on a 
self-consistently determined order parameter field, the situation would be even more com- 
plicated. As in the preceding section, the SN-action may also be subjected to a mean field 
approach. However, due to the imposed inhomogeneities in the order parameter, the sta- 
tionary phase configurations will in general no longer be spatially uniform. At first sight it 
may not be obvious how solutions to the spatially inhomogeneous stationary phase problem 
may be found. The correct strategy for this problem is again prescribed by the existence of 
the threefold hierarchy of energy scales, discussed above. Before going into details, let us 
give a brief account of the forthcoming construction steps: 

1. We will first employ the general ansatz 



an eigenvalue structure set by the matrix erf ® thereby automatically solving the 
saddle point problem with regard to the highest energy scale /1 (cf. the corresponding 
remarks made above). 
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2. In a second step we substitute the above ansatz into (|18D and derive a 'medium- 
energy' effective action that contains no energies higher than e and/or A. Thirdly 
we will perform a second stationary phase analysis thereby determining those field 
configurations (|30|) that extremise the medium energy action. 

3. By accounting for fluctuations around these configurations, we will finally be able to 
explore the low energy physics on scales u + . 

Beginning with the derivation of the 'medium energy action', we now formulate this program 
in more detail. 



C. Gradient expansion and 'medium energy action' 

In constructing the effective medium energy action, it is again crucial to exploit the 
existence of a scale separation in energy. Anticipating that the relevant field configurations 
T(r) fluctuate weakly as a function of r, we first borrow a parameterization of the kinetic 
energy operator that has previously been used in constructing the quasiclassical equations 
of superconductivity [01 



2m 



P 



-Aaf 

c 6 



a 



1 1 



1 



1 



— p + — p • q, 

2m m 



(31) 



where 



—id — 



-Aaf ® af 



The idea behind Eq. ([31]) is that the slowly fluctuating entities in the action, most notably 
Q, effectively do not vary on scales of the Fermi wavelength. Thus, it makes sense to 
decompose the momentum operator into two parts, p = p + q, where the 'fast' component, 
p, has eigenvalues of order of the Fermi momentum, pp, and so can be treated as a c-number 
with regard to slowly varying structures. The 'slow' component, q, accounts for both slow 
spatial variations and the magnetic field. For a more substantial discussion of ([31]) we refer 
to the original literature 0. 

We next substitute (130) and (BID into the action (HH) to obtain 



S[Q] 



-str r „ln 

2 ,p 



H - —p 2 + —af <g> of +V 1 + V 2 + T -1 [Vi + V 2 , T] 



2m 



2t 



(Go)- 1 



where 



V x = p • q, V 2 

m 



a 



P h 



UJ 



and str riP denotes a trace of internal indices, the 'fast' p's and the spatial coordinate. We 
next expand to lowest non- vanishing order in the 'slow' operators Vi. As will become clear 
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from the structure of the resulting series, the small parameters of the expansion are (l/L) 2 , 
for Vi, and er, Ar for V 2 - Here L denotes the typical scale at which the matrices T fluctuate. 
To lowest order we obtain 



S[Q) 



1 

— str r p 



G Q T- l [V 2 ,T) 



1 



G (^i + T-^T]) 



Note that there is no contribution at first order in V\. The reason is that V\ is linear in 
the vectorial fast momentum p, whilst Go is even in p. Thus, the trace over fast momenta 
annihilates this contribution. 

To prepare the tracing out of the fast momenta, we next formulate some useful identities 
describing the behaviour of the 'fast' Green function Gq. The following relations can be 
proved straightforwardly by explicitly performing the momentum integrations (cf. Eq. (|23|) ) 
and using some Pauli-matrix algebra. 



Go in its momentum representation may be written as (cf. Ref. 

1 



Go(p) 



E 

s=±l 



1 + saf 



< (To 



<ip) + sip 



and Eq.(|| 



The momentum trace over a single Green function becomes 

J2 G Q (p) = const. • 11 - mv n (Yo\)al h <g> a™. 



Further, if operators A and B vary slowly in space, then 



£str [Go(p)p-AG (p)p-B 



m 
T 



(Vol)2vrDz/ n ^str 



;i + saf )k ■ (1 - saf)B 



An application of these identities to the effective action above leads to 

S — >• S% + 5*2, 

where 



Si = —%- 



.-KUr, 



str 



Qaf® (A(r) + 6 + ^a : 



and, setting O = q + T 1 [q, T] — T 1 qT, 



S 2 = --7rDv n I str (1 + af ® a 3 ar )C(l - of <g> af)0 



By using the identity 



str 



P h 



-str 



[q,Qf 
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we obtain 



So 



7T 



str 



0Q8Q 



where 



d = d — i-Ala* 

c 



ph _ tr 



a 



3 • 



(32) 



represents the covariant derivative. Putting everything together we obtain our final result 
for the 'medium energy action' 



S[Q] 



7TZ/„ 



str 



D(t)Q) 2 + AiQaf ® f A + e + ^ai:' 



(33) 



Correlation functions are now obtained by substituting the action ([£J) into a functional 
integral over all Q-fields subject to the constraint Q 2 = II : 



S[Q] 



(34) 



where T>Q denotes the invariant measure on the manifold of matrices Q 2 = 1. For instance, 
the correlation function (13) takes the form (cf. Eq. (|I~7j)) 



(tr ph , r (£ r (e + u + /2)) tr ph>r (£ a (e - U+/2))) 



V 4 



VQ J str (g^a 3 bf ® af) J str (QE^af ® of) 



-S[Q] 



(35) 



In the limit A — > 0, the functional integral represents a superposition of two independent 
copies of normal metal u-models, one corresponding to the particle, the other to the hole 
sector. Due to the decoupling of these two components, the ph-structure becomes mean- 
ingless. For A ^ 0, the situation is more interesting. Given the spatially inhomogeneous 
structure of the action, an exact computation of correlation functions - in the sense of a 
complete integration over the nonlinear field manifold - is in general not feasible. Under 
these circumstances, the first and seemingly straightforward approach one might try is a per- 
turbative one. Yet, as usual with perturbative approaches in SN-physics, straightforward 
perturbation theory does not work here. 

To understand the origin of the difficulties let us introduce the parameterization 



Q 



of)e 



-XV 



(36) 



where [W, of 11 g) crf] + = 0. The parametrisation (|36| ) is frequently used in perturbative 
analyses of the a-model. In standard (N) applications of the a-model, its substitution into 
the action leads to a series 



S[Q] -> S[W] = S (2) [W] + S (4) [W] + S (G) [W] + 



(37) 



where S ,( ^ 2n ' ) [H / ] denotes the contribution of 2n-th order in W. The functional can then be 
evaluated by expanding perturbatively around the second order contribution exp(— S^fl^]) 
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and applying Wick's theorem. The resulting Taylor series converges rapidly due to the fact 
that contributions S"( 2n ) [W] are multiplied by large coupling constants g n 3> 1 (all of which 
are parametrically of the same order.) 

In the case A ^ the situation is more complicated. The point is that contributions 
S^[W], m being odd, arise from the perturbative expansion of the vertex ~ stT(Qcr^ h A). 
In particular, a non-vanishing contribution of first order in W emerges. The presence of 
this term invalidates perturbation theory. (This can be seen formally by means of a simple 
power counting argument: In the expansion of exp(— S^fVF]), each W is multiplied by a 
large coupling constant gi. On the other hand, the Wick contraction of two W's, gives 
a factor g^ 1 . Thus, the perturbative series expansion of exp(— [W]) diverges in the 
parameter g\j ' g 2 1.). 

To get some idea of how these problems can be overcome, it is helpful to understand 
the physical origin of the divergences arising in perturbation theory. To this end let us 
consider the disorder average of the ph-block of the Gorkov Green function, Q l2 (e), as a 
simple example of a quantity that strongly couples to the divergence of the W-perturbation 
series. When expressed in terms of the functional integral, the average (Q 12 (e)) takes the 
form 

«7' 12 (e;r,r)> ~ (str (Q(v)E^ ® E$ ® a?)) Q ~ (str(W(r)X)) w 

where (. . -)q/w stands for functional integration in the Q- respectively PF-representation of 
the theory and X is the fixed matrix, X = <g> E^ <2> erf. 

Suppose now, we intended to evaluate this functional expectation value perturbatively. 
To lowest order in W we would obtain 

(£ r ' 12 (e;r,r)) ~ ^str (W(r)X) str J dv'(W{r')K{v'))^ ~ J dr'K(r, r')A(r'), 

where (. . .)o stands for functional integration weighted by the quadratic action i>( 2 )[W] ~ 
/ str {W K~ l W). The kernel, K, governing S^jW] is the familiar diffusion pole K^ 1 ~ 
Dd 2 + ie. Thus we see that the first correction to the 'anomalous' Green function Q 12 (e) is 
proportional to the order parameter and - owing to the spatially long ranged behaviour of the 
diffusion pole - stretches far into the normal metal. Moreover, since the characteristic energy 
scale of the diffusion pole is max(e, E c ), we see that the correction is of 0(A/max(e, E c )), 
which, for sufficiently strong order parameter /coupling between N and S, exceeds unity. 
Remembering that the characteristic scale of the dimensionless quasiclassical Green-function 
is unity we have to conclude that the perturbation series resulting from a naive iy-expansion 
of the functional integral does not converge. To understand both the reason for this failure 
and the particular form of the first order correction, it is instructive to compare with the 
type of divergencies that appear within diagrammatic perturbation theory. In diagrammatic 
analyses, the correction to first order in A to the anomalous Green function is indeed given 
by a single diffusion mode. The real-space representation of this term has already appeared 
in fig. |](&), while the corresponding diagrammatic representation has also appeared as fig. |5|. 
Noting that this correction is only the first contribution to what becomes upon summation 
a full representation of the proximity effect, the origin of the problem becomes clear: By 
perturbatively expanding around af (8) erf r we have chosen the metallic limit of the Gorkov 
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Green function, Q a,r = iicf as a reference point. The superconductor, however, drives the 
adjacent normal metal region to a state that is far from conventionally metallic. In order 
to force a description of the system in terms of a perturbation theory around the metallic 
limit, we have to pay the price of an infinite perturbation series. Even worse, due to the 
effective spatial inhomogeneity of each perturbative contribution, arising from the space 
dependence of the diffuson, summation of the series becomes impossible. To summarize, 
the considerations above tell us that perturbative approaches based on spatially constant 
reference configurations are doomed to fail and that the origin of the problem is the spatially 
inhomogeneous manifestation of the proximity effect. 



VI. STATIONARY PHASE ANALYSIS 

Given what has been said at the end of the previous section, the correct strategy for 
overcoming the problems arising in perturbation theory becomes apparent: Prior to any 
perturbative attempts, it is preferable to seek a solution to the stationary phase equation 
5S[Q]/5Q = 0. Due to the spatial inhomogeneity of the problem, no uniform solutions 
Q = const, will be found. Once a solution Q has been obtained, both perturbative and 
non-perturbative evaluation schemes may be safely superimposed. The reason is that, by 
construction, no linear terms appear when the action is expanded around Q. 

We find it convenient to formulate the stationary phase analysis in a gauge where the 
phase dependence of the order parameters at the S/N-boundaries has been eliminated, at the 
expense of introducing a vector potential in the bulk N-region. To be specific, we perform 
the gauge transformation 



Q(r) 



exp 



ph 



tr 



- + (p(r) ) af ® ct 3 



Q(r) exp 



7T 



7T + ( f( r l 
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P h 
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tr 



(38) 



where, within the superconducting region, tp(r) is equal to the phases of the order parameter, 
as in Eq. (|To]), and in the normal region can be chosen arbitrarily. Inserting the gauge 
transformed field into the action we obtain 



S[Q] 



7TZA-, 



str 



D(dQ) 2 + AiQY ■ a ph 



(39) 



where 

Y = zA(r)e! + (e + w+<7?/2) e 3 , (40) 
and the vector potential entering the covariant derivative has been transformed by 



A 



-A - -t*p 



(41) 



Note that the right-hand side of Eq. (|4T| ) may be interpreted as —2m times the superfluid 
velocity. To find the stationary phase equation, we introduce a small variation 



Q 



'Qe 



Q + [SW, Q] 
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into the action and demand vanishing of the contribution at first order in SW. A straight- 
forward calculation then yields the equation 

Ddi{QdiQ) -i [Q, Y ■ <x ph ] = 0. (42) 

The first step to analysing the general set of solutions of this equation is to specify that the 
solution is as simple as possible, i.e. as diagonal as possible. Noting that the equation is 
diagonal in ar-,tr- and bf-space (in bf it is even trivial), we see that, as in the case of a bulk 
superconductor, a sufficiently general ansatz reads 



bf,ar,tD 



(43) 



where q± refers to the retarded/advanced blocks defined in Eq. ( |21~D and q± 2 are vectors 
of complex numbers (i.e. structureless in bf-space). The restriction of the saddle point 
equation to the blocks q± = q± ■ a ph now reads 



where 



Ddi(q±diq±) 



a 



ph 



0. 



(44) 



zA(r)ei + (e±oo + /2)e 3 . 

Comparing ([TJ) with (|]) and identifying q± with g^, we identify the stationary phase equa- 
tion of the nonlinear cx-model as the Usadel equation. One consequence is that we are 
immediately able to write down the boundary conditions at a (perfectly transmitting) SN 
interface, by direct analogy with the Kuprianov and Lukichev relations, eqns. ( |10a| ) and 
(ED: 



<yq±d±q± 



q± 



o~q±d±q± 



q± 



(45a) 
(45b) 



the former of which implies current conservation at the interface. Here d± is the normal 
derivative across the (planar) boundary, and x± denotes a space point infinitesimally to the 
left/right of a boundary point x. Note that the normal state conductivities, cx, in the left 
and right region may differ. 

In passing we note that, although we have stated above and will use further the relations 
for a perfectly transmitting interface, in general we need not keep to such a restriction within 
this formalism. For instance, we could have equally well employed the following conditions 
in the limit of small transparency, again by analogy with eqns. ( |10a| ) and (|10b|) : 



vq±d±q± 



o-q±d±q± 
Gt r 



q± 



(46a) 
(46b) 



where Gt is the tunnel conductance of the junction, as given by Eq. ( |A5|) . By modelling the 
tunnel barriers microscopically (as was done, e.g. in Ref. |5T|) these boundary conditions 
can in fact be rederived within the cr-model formalism. 

The coincidence of the stationary phase equation of the cx-model with the Usadel equa- 
tion, which, as mentioned in the introduction, was first observed by Muzykantskii and 



Khmelnitskii [23] in a different context, has fundamental consequences for all that follows: 
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• Already on the level of the saddle point equation, the cx-model contains all the in- 
formation that is otherwise obtained quasiclassically. In particular, solutions of the 
equation can in most cases of interest be imported from the extensive literature on 
Usadel equations for SN-systems. 

• The facts that a) the solutions of the Usadel equation for the retarded and the advanced 
Green function are different and b) the uj + = action is isotropic in ar-space, imply 
that we encounter a situation of spontaneous symmetry breaking: The mean field does 
not share the symmetries of the action and a Goldstone mode, operating in ar-space, 
will appear. 

These two observations suffice to dictate the further strategy: One first has to solve or 
import a solution of Eq. (^). Then the solution Q = diag(g + ,g_) is substituted back 
into the action and fluctuations around the block diagonal solution are introduced via, 
Q — > TQT^ 1 . In analysing fluctuations, the main emphasis will be on exploring the role of 
the Goldstone mode. However, before we proceed to the actual formulation of this program, 
it is worthwhile to stay for a moment at the mean field level and to acquire some familiarity 
with the Usadel equation and the structure of its solution. 

We first note that the different 'sectors', q± 2 , of the solution vector are not independent 
but rather connected to each other via symmetry relations. 

1. The general relation (cf. Eq. (g)) 

g A (v,r')=af(g R (v',v)yaf (47) 

implies 

q_ = diag(l,l,-l)(q + )*. (48) 

2. Taking the transpose of the Usadel equation in the tr-sector 1, we obtain 

q 2 = diag(l,-l,l)q 1 . 

As for the spatial behaviour of the solution, some remarks may be made in general. Deep in 
a superconducting region, the large cxf 11 - component of Y enforces an approximate equality of 
q ~ e.\. Conversely, deep in a normal metal, q will be aligned with The Usadel equation 
describes a smooth interpolation between these two limits, where the gradient term inhibits 
strong spatial fluctuations (cf. fig. ||). 
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FIG. 6. Schematic plot of the retarded component (real part) of the saddle-point solution, 
q_l_, for an SNS-j unction with coincident phases of the order parameters. 

In order to say more about the spatial structure of the solution to the Usadel equation, 
we have to restrict the discussion to specific examples. Here we will consider two simple 
prototype systems, representative of the wide classes of systems with a) quasi infinite, and 
b) compact normal metal region. Since we consider a quasi- ID geometry in each case, we 
denote by x the position variable perpendicular to the interface. 

Infinite SN-junction: Consider the model system shown in fig. [7|. The normal metal and 
superconductor regions occupy x > and x < respectively, so that the gap function is 
modeled by A(x) = |A|G(— x). The system is quasi one-dimensional in the sense that its 
constant width is comparable with the elastic mean free path (i.e. there are many conducting 
channels but no diffusive motion in the transverse direction.) We assume that no external 
magnetic field is present. Furthermore, since there is only one superconducting terminal, an 
elimination of the phase of the order parameter does not induce a gauge potential and we 
can globally set A = 0. 




FIG. 7. The geometry of the infinite, quasi-lD SN junction. 



The analytic solution of the corresponding Usadel equation is reviewed in Appendix |U 1| . 
Due to the global absence of a vector potential, the spatial rotation of the vector q takes 
place in the 1-3-plane only. Hence, it can be parameterized as (cf. the analogous form for a 
bulk superconductor, Eq.(p6|)) 

q(x) = (sin9(x), 0, cos9(x)). 

In fig. H] we have plotted the curve in the complex plane that is traced out by 9(x) upon 
variation of x. 

IrrnH 



71/2 Re-d 
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FIG. 8. Dotted line: the trajectory of the retarded component of 9 as a function of position for 
an infinite SN junction with coincident phases and e/A = 0.6 and 7 = 0.1. The trajectory starts 
at x = and approaches the origin as x — > 00. Solid line: The same for an SNS system, with 
L/£ = 5, starting at x = —L/2 and ending at x = L/2. The trajectory reverses direction at the 
centre of the N region, x = 0. 

Fig. y shows the local DoS z/(x) = z/ n Recos6 l (a;) (cf. Eq. (|28|)) obtained for a particular 
value of the material parameter 7 = v n \fD n / (z/ s y^DI), where D n)S are the diffusion constants 
in the N,S-region. Note that the deeper one proceeds into the S-region, the more the DoS 
approaches the characteristic BCS-form. Due to the proximity effect, the structure of the 
subgap (e < |A|) DoS in the N-region remains non-trivial. Only in the asymptotic limit 
x — > —00, the region of suppressed DoS shrinks to zero and 'normal' behaviour is restored. 
More precisely, substantial alteration of the DoS induced by the proximity effect is restricted 
spatially to a region of several diffusion lengths from the interface into the normal metal, or 
several (dirty) superconducting coherence lengths into the superconductor. 




FIG. 9. Local DoS of the infinite SN junction as a function of both energy and position for 
7 = 0.1. 

Fig. shows how variation in 7 affects the DoS. In particular, we may take the 'rigid' 
limit 7 — > 0, for which the bulk superconducting value of the Usadel angle is imposed at the 
interface, and retain a non-trivial structure in the spectrum. 
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FIG. 10. The effect of 7 on the energy dependence of the DoS, at x = 1.5^. 



For a more comprehensive discussion of the local DoS of the system above we refer the 
reader to e.g. Ref. [fT(i|l . 

SNS-junction: As an example representative for the class of SN-systems with compact 
metallic region we next discuss the quasi one- dimensional SNS-junction displayed in fig. ||. 
As opposed to the SN-system, the physics of the SNS-system does depend on the phases of 
the order parameters in the superconducting terminals. For this reason SNS-systems with 
an in-built possibility to tune the phases of the superconductors are sometimes referred to as 
Andreev interferometers. In Appendix |C 1| the quantitative solution of the Usadel equation 
is discussed explicitly for the comparatively simple case Aip = ipx — <p 2 — 0. 

An analytic solution to the general case, A(p ^ 0, is also possible, although cumber- 
some. In the limits of a very short (L <C L e ) and very long (L ^> L e ) junction, separate 
approximation schemes to the general solution may be employed: the short junction case 
has been treated by Kulik and Omelyanchuk |36| (for L <C £, L e ) and Likharev |]30|j32 



(for ( C L < L e ), while the long junction case (L ^> £, L e ) has been treated by Zaikin 



and Zarkov [31j33|. At this point we restrict ourselves to a discussion of a few qualitative 
characteristics of the solutions obtained in various regimes. 

The most conspicuous feature of the system is the appearance of a minigap, E g . As 
mentioned in section [0], the precise form of the gap depends on the phase difference Aip. 
Its maximum width is taken for 5ip = 0, while this width tends to A and E c for short and 
long (as defined above) junctions, respectively. With growing phase difference, E g shrinks 
until, at Aip = 7r, it closes altogether [[J^J (apart from a 'microgap' of the width of the 
single particle level spacing |10[ - the latter is induced by the general phenomenon of level 
repulsion in disordered metals). 

We take as illustrative the case of an SNS-junction with coincident phases, and of arbi- 
trary width. In this case, the solution for q again lies in the 1-3-plane and is parametrized, 
as before, by the angle cos 9 — q • £3. Fig. || shows the trajectory of the retarded 9 in the 
complex plane as a function of position, and fig. |TT| shows the local DoS as a function of 
both energy and position. Note, as compared to the infinite SN case (fig. |9p, the appearance 
of the minigap, E g , which is below the superconducting gap and independent of position. 
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FIG. 11. The local DoS for the SNS-junction as a function of both energy and position, for 
L/£ = 5, 7 = 0.1 and if = 0. 



Figs. and [13| show how variation of the material parameter 7 affects the local DoS. 
Further, fig. [13] shows how the local DoS at the SN-interface decreases as 7 is reduced (e.g. 
as the ratio of the disorder concentrations in the S- and the N-region is lowered.) 



> 
> 



2.0 



1.0 



0.0 



Y=0 




_ y 1 




- Y 5 






/ ~ ~\ ~~ — 






1 , 
1 1 
/ , 

1 

1 1 
1 1 
1 1 

1 ; j 







0.0 



0.5 



1.0 



e/A 



FIG. 12. The effect of 7 on the energy dependence of the local DoS at the centre of the junction. 
Here L/£ = 3. 
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FIG. 13. The effect of 7 on the position dependence of the DoS for the SNS-junction, with 
e/A = 0.5 and L/£ = 3. 



VII. FLUCTUATIONS 

Having discussed the mean field content of the theory, we now turn to the central issue of 
this paper, that of mesoscopic fluctuations. The generality of the field theoretic machinery we 
have been discussing is such that it may be employed to analyse the majority of fluctuating 
observables in SN-systems. Our strategy of extending the conventional cr-model formalism 
allows us to take full advantage of the versatility of an approach that has been greatly 
developed in the study of fluctuation phenomena in purely N-mesoscopic systems. Yet in 
order not to diversify too much, we focus here on the discussion of the specific example 
of spectral fluctuations. In fact, the status of spectral fluctuations is slightly higher than 
that of an ordinary example, since a vast number of fluctuation phenomena are directly or 
indirectly related to fluctuations in the single particle spectrum (see Ref. p4|,^8| for review). 



A specific issue is the nature of fluctuations in the local DoS around the mean values 
(displayed in figs. |9| and [lTJ) . To characterize these fluctuations quantitatively, we employ 
the correlation function (0). Representing the correlation function in terms of Green func- 
tions (cf. Eq. (0)), we see that we need to compute functional expectation values of the 
type specified in Eq. ([35]). An evaluation of Eq. (|35| ) on the mean field level discussed pre- 
viously leads to a vanishing fluctuation component, (G r Q & ) — (G T )(Q B '), due to the fact that 
the mean field configuration, Q, is diagonal in the ar-indices. In other words, there is no 
connection between advanced and retarded components and the functional evaluation of the 
product of Green functions equals the product of the averages. Consequently the physics of 
mesoscopic fluctuations is contained entirely in fluctuations around the block diagonal mean 
field configuration: Q — > TQT^ 1 . At first sight the analysis of the fluctuation degrees of 
freedom, T - after all a 16-dimensional matrix - seems to be a complicated task. Fortunately 
the totality of fluctuations may be organized into three separate types, each with its own 
physical significance. Such a classification scheme leads to a substantial simplification of the 
analysis. To be specific, we distinguish between fluctuation matrices, T, that 

a) are diagonal in the space of advanced and retarded components. These types of fluc- 
tuations still do not give rise to correlations between retarded and advanced Green 
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functions. Nonetheless they are of physical significance: Quantum corrections to the 
Usadel solution are described by fluctuations of this type. We elaborate on these effects 
in section |V11 A| . 

b) are non-diagonal in ar-space but are proportional to unity in ph-space, [T, o~ ph ] = 0. 
Fluctuations of this type do induce correlations between different Green functions and 
thereby mesoscopic fluctuations. They will be discussed in detail in section |V11B. 



c) fulfill neither of the conditions a) and b). Whereas the physical significance of these 
fluctuations is less clearly identifiable than the one of the a- or b-type fluctuations, 
they are nevertheless of importance. The reason is that the c-type fluctuations tend 
to destroy at sufficiently high energies the correlations that derive from mesoscopic 
b-type fluctuations. 

The three fluctuation types are summarized in the following table: 



type 


structure in ar space 


ph space 


lead to 


a 
b 
c 


diagonal 
non-diagonal 
non-diagonal 


non-diagonal 

oc ll ph 
non-diagonal 


corrections to Usadel Green function 
mesoscopic fluctuations 
destruction of b-type fluctuations 



Besides the criteria a)-c), further restrictions to be imposed on the fluctuation matrices 
follow from two fundamental symmetries of the model: Firstly, general convergence criteria 
|64]| enforce the condition 



T f = r]T- V\ 



(49) 



where 



n 



E 



bf 
11 



a 



P h 



On -f- tLl, 



bf 
22- 



Secondly, the tr-space structure of the model implies |64| 

T T = tT- 1 t-\ 

where 

iaf + E\{ ® a? . 



(50) 



^22 



For future reference we note that it is often convenient to represent both the matrices T and 
the above symmetries in terms of the generators of the fluctuations: 



T = exp(W), 

where the generators W are subject to the constraints 
and 



(51) 



(52) 
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W 1 = -tWt-\ (53) 

The above scheme generally classifies the various types of fluctuation corrections to the 
quasiclassical picture of dirty superconductivity. Beginning with the a-modes, we now turn 
to a more comprehensive discussion of these fluctuations. Although we do not provide 
by any means a comprehensive survey of the full diversity of effects, we will find even in 
simple geometries of SNS-structures a significant range of fluctuation types and associated 
phenomena. In order to prevent the subdivision of these phenomena into excessively many 
classes, we assume throughout that the spatial extent of the N region is sufficiently large 
that E c <C A. In this case, the minigap of the SNS-junction with zero phase difference, 
A(p = 0, lies at energy E c . The significance of this restriction will be discussed further 



below, in section VIII. 



A. a- Type Fluctuations: Quantum Corrections to the Quasiclassical Theory 

In this section, fluctuations of type a) around the Usadel saddle point will be considered. 
After specifying the general structure of these fluctuations, we will exemplify their effect by 
discussing the quantum corrections to the quasiclassical picture of single particle properties 
of SNS-structures. Other types of SN-structures can straightforwardly be subjected to an 
analogous analysis. 

Fluctuations of type a) are diagonal in ar-space. Since the saddle point Q is also ar- 
diagonal, the full effect of the a-type fluctuations may be studied by considering just one 
of the ar-sectors of the model action fl33|). For example, we may consider the retarded 
sector describing the behaviour of a single, disorder-averaged retarded Green function. The 
restricted action is given by 

S re t[Q n ] = / str [D(dQ ll Y + 4*Q n (z|A|a ph + ea 3 ph )] , (54) 

where the eight dimensional matrix field Q u denotes the rr-block of Q and the parameter uj 
(which is meaningless for a single Green function) has been dropped. To keep the notation 
simple, we will henceforth (until the end of the section) denote Q 11 again by Q. 

Following the general philosophy of our classification, we organize the field Q into a 
saddle point contribution Q (which is given by the solution of the retarded Usadel equation) 
and a-type fluctuations around it: 

Q = RQ f R~\ Q f = TafT~\ (55) 

where R represents the inhomogeneous rotation parametrising the saddle-point, 

Q = q. (T ph = Rag* BT 1 , (56) 

and the rotation matrices T generate the a-fluctuations. More precisely, T G G/H, where 
G is the group of eight dimensional matrices subject to the constraints (see Eq. fl49|)), 

Tt = Va T- l V :\ Va = < ® of + EH (57) 
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and Eq. (|50|). The subgroup HcGis denned through H = {h G G\[h,af] = O}0. Note 
that these symmetry relations imply that the matrices T span Efetov's eight dimensional 
coset space of orthogonal symmetry. In other words, save for the presence of the order 
parameter (and the different coupling of magnetic fields), the action (|54j) is identical to that 
of a conventional cr-model for an advanced and a retarded normal metal Green function. 

In order to understand more fully the effect of the fluctuation matrices T, we first have to 
analyze their commutation behaviour with the different contributions to the action (|54|), i.e. 
the magnetic field, the energy term and the order parameter. Surprisingly, it turns out that 
there is a subset of T's which not only commute through the order parameter (no matter 
what its phase) but also are insensitive to magnetic fields. For reasons that will become 
clear below, we will call these matrices the C-modes. In the limit e — > 0, the C-modes 
become completely massless. This implies that for low energies these modes need a special, 
or non-perturbative, treatment. 

1. The C-Modes: Non-Perturbative Corrections to Quasiclassical Green Functions 

Amongst the set of a-type matrices, T, consider the subset T c = exp(Wc), subject to 
the additional constraints, 

1. [cf, l^c] = (order parameter 1-component commutes through), and 

2. [af 11 (g> af, Wc) = (no coupling to magnetic field). 

In combination with eqns. fl52|) and (|53|), this gives altogether four constraints and the non- 
trivial statement is that a set of generators Wc obeying all of them actually exists. These are 



Note that there is some freedom in parameterizing fluctuating field configurations. For example, 
as an alternative to Eq. (|55|), one might let the fluctuation matrices T act from the 'outside', 

Q = TQT- 1 . 

As we will see shortly, there are situations where this parameterization is advantageous. In general, 
however, it creates unwanted problems. To see this, interpret the T's as 'rotations' acting on the 
unit vector q appearing in Q = q • cr ph . Clearly, there are rotations that are ineffective (namely 
those that rotate q around itself) and should be excluded from the configuration space of the 
T's. In practice, however, it is difficult to disentangle these rotations from the relevant ones. 
For example, parameterizing the matrices T in terms of some kind of spatially fixed coordinate 
systems, T = T(9±, . . .), where 9i are rotation angles around certain fixed axes, one finds that the 
effective action S[9i, . . .] contains unphysical divergencies. These are due to the fact that some 
'directions' in the parameter space spanned by the 0's correspond to ineffective rotations, thereby 
being energetically costless. The way to remove these spurious degrees of freedom is to introduce 
a 'moving' coordinate frame which, by construction, only parameterizes rotations around axes 
perpendicular to q. This is exactly what the representation Q55| ) achieves. Due to T £ G/H, the 
T's contain only degrees of freedom that effectively modify the matrix a^ h . 
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the C-modes. Before turning to the actual construction of these modes, let us qualitatively 
discuss some of their general properties. 

First note that the conditions 1. and 2. above imply that the C-modes further commute 
with the 2-component of the order parameter. This follows from the observation that the 
phase twist needed to interpolate between the 1- and the 2-component of the order parameter 
is equivalent to the appearance of a magnetic gauge field which (see 2.) is invisible to the 
Wc's. Being insensitive to both magnetic fields and order parameters with arbitrary phase 
positioning, the C-modes merely couple to the gradients and the energy term in the action 
(f>3|). In the limit of small e, they become completely massless. More precisely, for a spatially 
constant T® 6 Gc, where Gc is the subgroup of G fulfilling the extra constraints 1. and 2., 

S[T°Q(T°rV=°S[Q]. 

Physically, the C-modes represent modes of quantum interference in SN-systems which 
survive magnetic fields. Within a different formalism, these modes have for the first time 
been noticed in Ref . |T7| . Subsequently various physical phenomena caused by their presence 



have been discussed in the literature: 

• In the initial paper, Ref. |T7| mentioned above, it was observed that in SN-systems weak 



localization corrections to the conductance may survive magnetic fields. The quantum 
interference process responsible for that effect is associated with the C-modes. 

In Ref. [|1(J, an SN-system subject to a magnetic field, but not exhibiting a minigap 
(due to suppression of the proximity effect by the magnetic field), was considered. In 
this case the single-particle DoS vanishes at the chemical potential on a scale set by 
the mean level spacing. The existence of this 'micro-gap' is also an effect caused by 
the C-modes. 



Feigl'man and Skvortsov |6l] discuss the effect of C-modes on the transport behaviour 
of vortices in moderately clean type II-superconductors. 

The above phenomena relate to mean single particle properties. For the sake of com- 
pleteness, we here mention some manifestations of C-mode fluctuations in two-particle 
properties: The level statistics of SN-quantum dots in a magnetic field falls into a 
symmetry class that is different from any of the standard Wigner-Dyson classes. Re- 
ferring to a classification scheme due to Cartan, the SN/magnetic field symmetry class 



has been termed 'class C pOj] . As in standard mesosocopic systems, these level fluc- 
tuations can also be associated with channels of microscopic quantum interference. 
Whereas Wigner-Dyson fluctuations in diffusive N-systems are caused by 'diffuson' 
and 'Cooperon' modes, the class C fluctuations are connected to the modes specified 
above, and hence the name 'C-modes'. A u-model formulation of the C-mode spectral 
statistics of random matrix ensembles was presented in Ref. [19| . 



C-type level statistics in vortices has recently been microscopically derived by 
Skvortsov et al. []62 |. Thermal transport carried by C-modes through the core of 



superconductor vortices is considered in a recent paper by Bundschuh et al. p3 



Being effective already on the level of single particle properties, the C-modes must originate 
from interference processes between particles and holes. However, they cannot be identical 
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with the modes displayed in figure p] (b), since the latter are field sensitive. A typical type-C 
path configuration is displayed in figure [14|. In the analysis below we will derive quantitative 
expressions for processes of this type. 




FIG. 14. Semiclassical illustration of an interfering path configuration contributing to the 
C-mode corrections to the DoS. Compare the relative orientation of the arrows with those ap- 
pearing in figure [l] (b). 

After these general remarks, we next turn to the analysis of the C-mode contribution 
within the cx-model formalism. Specifically, we will discuss the effect of these modes on the 
mean density of states of the SNS-geometry displayed in figure 0. 

Before turning to the core of the discussion, let us make a technical remark which will 
have some impact on the organization of the remainder of the section: The invariance of 
the action under spatially constant C-transformations , S[T®Q (Tq) ] e = S[Q] suggests an 
interpretation of these modes as a global symmetry of the action. In particular, it is more 
natural to let them act on Q from the outside: Q = TcQT^ 1 (compare with the inside 
representation of general a-fluctuations, Eq. (|55|) , and the footnote on p.f|0]). Of course 
it is possible to forcefully contrive an inside parameterization for the C-fluctuations, via 
TcQTq 1 = RTc&zTc 1 Rr l , where the unitarily transformed T = R~ l T c R- In practice, how- 
ever, this representation is inconvenient and, more seriously, makes it difficult to separate 
the C-mode fluctuations from general fluctuations around the saddle point. These consider- 
ations imply that, in general, it is difficult to treat the C-modes and the rest of the a-type 
fluctuations simultaneously. Physically, however, these problems are of little significance: 
Below it will be shown that the minimum price in energy associated with a non-C-fluctuation 
is of 0(E C ). This implies that two regimes with qualitatively different fluctuation behaviour 
exist: 

• Low energies, e<£ c , where the C-modes are relevant, whereas the other a-fluctuations 
can safely be ignored, and 

• High energies e > E c , where all fluctuations have a comparable action of 0(E C ) and it 
is pointless to carefully distinguish between the different types (C or other). 

Below we will discuss these cases separately. Although our analysis is not applicable to 
the crossover regime of intermediate energies, we do not expect qualitatively remarkable 
phenomena to arise there. 

Low energies, e C £ c : To quantitatively analyse the fluctuation physics in this regime, 
we need first to derive an effective action for the C-modes. Fluctuations other than C are 
ignored. 
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To this end, we use that [To, crf h ] = [T c , cr^ 11 ] = and represent the Q-Held as 

Q = qsQc, Qc = T c afT c \ (58) 
Substituting this parameterization into Eq. (|54"D, we obtain the desired action 



Sc[Q 



7TZ/„ 



C 



str 



D q 2 3 (dQcY + 4tq 3 eQ c a p 3 



P h 



(59) 



Among the general set of C-type fields Qc , there is a particular mode Q c = T c a 3 h (T c )~ l 
which not only has C-symmetries but also is spatially constant. Substituting this 'zero-mode' 
into (|59|) , we obtain the action 



SclQc] 



is 

— str 
2 



where 



(60) 



(61) 



Note that deep in N, where q 3 = 1, s = 7re/d, where d is the mean level spacing, coincides with 
the standard parameter s [Q commonly employed in the literature on spectral correlations 
in metals. 

In order to give the above zero mode action some physical significance, it has to be shown 
that it is energetically gapped against the action of the higher (spatially fluctuating) field 
configurations of C-symmetry. The spectrum of fluctuating C-modes can be determined at 
various levels of accuracy. For our purposes, it suffices a) to demonstrate that a ground 
state gap exists and b) to coarsely estimate its magnitude. To do so, we first note that the 
mode spectrum is essentially determined by the gradient operator appearing in Eq. 
Integrating by parts, the latter can be rewritten as 



?3 Str 



QDqfd{qld)Q 



(62) 



The rationale behind this reformulation is that Dq 3 2 d(q 2 d) may be regarded as a differential 
operator which is Hermitian with respect to the scalar product (/, g) = / q\jg- Being 
Hermitian with compact support (taking q 3 — > in S), the spectrum of the differential 
operator is discrete. To estimate the spacing AE between the zero eigenvalue of the spatially 
constant eigenmode and the first excited eigenvalue, we use the fact that the range of 
support of the operator is set by L, the extension of N. Standard reasoning for the eigenvalue 
structure of one-dimensional Hermitian differential operators with compact support then 
leads to the estimate AE ~ D/L 2 = E c . Notice that the actual spatial structure of the 
excited eigenf unctions may be complicated. For example, unlike with standard applications 
of cr-models to N-systems, typical eigenmodes of the action obey neither Neumann nor 
Dirichlet boundary conditions, but rather exhibit more complicated edge behaviour which, 
in principle, may be derived from Eq. ( |62|) once ^3 is known. For our purposes, entering this 
discussion will turn out to be unnecessary. 
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The considerations above show that for energies e <C E c , the action is governed by 
the spatially constant C-modef^- It is reasonable to ask for which physical applications 
such low energies may be expected to play a role. For zero phase difference between the 
superconductor terminals, the minigap is of 0(E C ) and, with regard to the DoS, the C-mode 
fluctuations are expected to be of little importance. To actually make visible the impact of 
these fluctuations on the DoS, we concentrate here on the case of a junction close to 7r-phase 
difference, i.e. a junction where the gap is nearly but not completely closed: E g <C E c . Note 
that previous studies of the C-type fluctuations having concentrated on the limiting case 
where the proximity effect is totally suppressed (E g — > 0). It is worth remarking that, while 
we have limited here our comments to issues surrounding the DoS, subgap properties of SNS- 
junctions, such as the Josephson coupling, may well be affected by the C-fluctuations, even 
under the broader conditions of a fully-established proximity effect. However, the analysis 
of these phenomena lies beyond the scope of this paper. 

Focusing on the range of applicability of the action ([Jl]), we notice that for very small 
energy parameters, so that s = 0(1), the zero mode must be treated in a non-perturbative 
manner. Very much as in the study of N cx-models at low energy scales, the fluctuations be- 
come unbounded as s — > 0. Rather than perturbatively expanding in terms of the generators 
Wc, it then becomes necessary to integrate over the entire manifold of matrices Tq- Non- 
perturbative analyses of this type have previously been applied to the study of a random 
matrix ensemble |19| of non-proximity effect SN-structures and of normal core excitations 
of vortices in superconductors |)2]| . Here we discuss how the C-fluctuations affect the low 
energy behaviour of the Gorkov Green function in SNS-junctions. 

To this end, we first need to specify a global parameterization of the matrices Tc- It is 
a straightforward (if lengthy) matter to show that a general 8 x 8-matrix Tq = exp(Wc), 
subject to both the general constraints Q52p and fl53p and 1. and 2., can be parameterized 
as 

(63) 




Where A and /i are Grassmann variables. For the (invariant) measure associated to the 
integration over the matrix Qc we obtain (cf. the analogous but more difficult calculations 
of integration measures in ||64|| ) 



,^ / x f 2n dy r d6 sin 9 . , , , 

dQc(. ■■)= 7T o • Vnio\ / dXd ^ ■ ' ' ' 
Jo 2ir Jo 2 sin (y/2) 



12 A closer analysis, similar to that presented in section VII Q with regard to the destruction of 
correlations by c-type fluctuations, shows that the range of stability of the C-zero mode is in fact 
limited by e ~ dy/g <C E c . 
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Substituting the parametrisation (|63| ) into the action (|60[) we obtain 



S° C [Q 



c\ 



-2is(cos9 - 1) 



(64) 



As an example, we apply this action to a calculation of the C-mode corrections to the local 
DoS. Substituting the zero mode integration over Qc for the functional expectation value 
in Eq. (|28| ) and performing the (trivial) integrations over Grassmann variables and y, we 
obtain 

v = u n Re |g 3 (l - ^ J* d6 sm6e 2lS{cose - 1) 



u n Re <q 3 1 



1 - e 



-4is* 



4is 



_ f / sin(4s) 1 - cos(4s) 
v n Re g 3 1 ^ + 



4,5 



4zs 



The last line tells us that for small energies, s — > 0, the DoS always (i.e. including the 
case of a 7r-junction) vanishes on a scale set by the mean level spacing. This is the DoS 
'micro-gap' that has been discussed previously in Refs. [|l^,|l^,^l|,|^] . Moreover, for general 
s, possessing real and imaginary components, the DoS is not only determined by Re q 3 
but also by the imaginary component, Im q 3 , of the Usadel solution. Finally, s contains 
the Usadel solution in an integrated form, that is, the C mode introduces some non-local 
influence on the local DoS by the Usadel solution at different points of the system. The 
corrections vanish algebraically as s -1 . For Im s > 1 the (oscillatory) factors containing the 
exponentiated parameter s can be neglected and we obtain the simplified result 



v 



4 s J 



(65) 



Note that the definition (61) implies that for the globalDoS (~ / v) the algebraic corrections 
vanish. For a diagrammatic interpretation of this correction to the Usadel DoS, cf. Ref. [ITU . 

For Aip = 71 — 0.0025 and g = 5, the quasiclassical DoS and the corrections to it are 
displayed in figure 15j. We display here the DoS in the vicinity of the minigap, at an 
energy regime similar to that considered by Zhou et al. f7I |. For this set of parameters, the 
quasiclassical DoS displays a very strong peak, which dies down for (e — E g )/E g = 0(1). 
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FIG. 15. (a) Quasiclassical DoS and (b) the C-mode correction for Aip = tt — 0.0025 and g = 5, 
at the centre of the junction and in the vicinity of the minigap. The dotted line represents the 
non-oscillatory part of the correction, as given by Eq. (|65|). 

Large energies, e > E c : For large energies, the isolated 0-mode action is no longer of 
significance. The energy of all modes, C- or other, is larger than E c . In particular, spatially 
fluctuating configurations (with a 'kinetic' energy cost of 0(D/L 2 = E c )) need to be taken 
into account, too. Since the energy associated with all these fluctuations is parametrically 
of the same order, separating the C-modes from the rest becomes pointless. In the next 
section we discuss the corrections to quasiclassics at energies larger or comparable to E c 
arising from a perturbative treatment of all a-type fluctuations. 



2. Perturbative Corrections to the Quasiclassical Green Function 



In the following we consider the impact of a-type fluctuations on the 'high' energy (e 
E c ) behaviour of the average Gorkov Green function. To keep our discussion simple, we limit 
consideration in this section to a SNS geometry where time-reversal symmetry is maintained 
(i.e. where the phase difference between the two superconductors, Aip, is zero.) The violation 
of time-reversal symmetry through the variation of the phase of the order parameter across 
the junction will not change our discussion qualitatively. 

Specifically, the questions we are going to address are 

• Do quantum corrections lead to the suppression of the minigap in the normal region? 

• If so, does the minigap edge remain sharp, or are states introduced at all energy scales 
below the gap? 

As in the previous sections, the theory developed in this section may also be straightforwardly 
generalized to other types of geometries and observables. 

For large energies e ~ E c , it is convenient to parameterize the totality of a-type fluctua- 
tions as in Eq. fl55p. Due to the comparatively large energy cost associated with fluctuations 
around the Usadel saddle point, it is sufficient to expand the action to low orders in terms 
of the generators W of the fluctuation matrices T. Substituting the parameterization ( |55|) 
into the action (0), we obtain 

S[Q] = -^p / str [DidQff + AieR^afRQf] , (66) 

where d = d + [R~ l dR, ]. Note that the above action does not contain the superconduct- 
ing order parameter. This is accomplished by demanding that the non-C-type fluctuations 
obey Dirichlet boundary conditions at the NS-interface. As a consequence these modes are 
spatially varying with a minimum fluctuation energy of 0(E C ), which justifies their pertur- 
bative treatment. As for the C-modes, these do not couple to the order-parameter anyway. 
In principle, their treatment is difficult because, as mentioned above, they a) fulfill mixed 
boundary conditions different from Dirichlet or Neumann and b) are difficult to separate 
from the complementary set of a-type fluctuations. However, for large energies we believe 
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these complications to be physically irrelevant: The spectrum of all fluctuating modes is 
discrete with a typical spacing of 0(E C ). For energies e comparable with E c , all modes need 
to be summed over. Under these conditions, the detailed structure of boundary conditions 
and/or eigenvalues of individual modes becomes largely inessential; What matter are the 
global features of the energy spectrum associated with the fluctuations, most importantly, 
the typical mean energy spacing. For this reason, we feel justified in ignoring the different 
boundary behaviour of the C-modes and to globally impose Dirichlet boundary conditions 
(thereby correctly modelling the typical spacing between consecutive eigenmodes). We be- 
lieve that this simplification does not lead to qualitative errors. 

To obtain the perturbative expansion of the action, we employ the exponential param- 
eterization T = exp(TU) and expand the generators W in terms of ph-Pauli matrices, 
W = Wl of + w 2 af (so that [W,af}+ = 0). For zero phase difference, the Usadel so- 
lution encoded in the rotation matrices R can be parameterized in terms of a single angle 6 
(cf. Eq. (|CID). The rotation matrices R mediating between af 1 and the Usadel saddle point 
then take the simple form 



R 



exp 



.0 P h 
2 2 



(dR)R- 



-atd9. 



(67) 



Substituting these expressions into the action and expanding up to second order in w iy it is 
a straightforward matter to show that 



S[W] = niy n f str [D ((dwi) 2 + (dw 2 f - (dOfw 2 ) - 2itcosd (w 2 + w 2 )] + 0(W 3 ). (68) 



Here and in the following, the notation 'stro' represents a supertrace over all degrees of 
freedom except for the ph-components, which have been traced over. The absence of terms 
at first order in W is assured by the expansion around the saddle-point configuration of Q. 

To eliminate the term in ~ (<96>) 2 , we make use of the fact that the Usadel equation ( |C2| ) 
possesses the first integral, 



D{d0) 2 - 42e(cosfl - cos 6(0)) = 0, 



(69) 



where we have used the fact that in the middle of the junction, 86 = 0. Substituting this 
result into Eq. flBlf), we obtain 



S[W] = 7rz/ n / str DUdwi) 2 + (dw 2 y)+ 2iewf (2 cos 6(0) -3 cos 6)- 2iew 2 i cos 6 . (70) 



To compute corrections to the DoS, we substitute the exponential parameterization into the 
functional representation (ESI) to find 



u( x ) = ^Re (str (ifcr* <g> afR^Q 



z/ n Re cos 6 
UnRe cos 6 



l + i(str (af(w\ + wl 



w 



nv n L 



(nifoz) - H-2(x,x)) 



(71) 
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where x is the coordinate along the junction and L± is spatial extent of the N region in all 
other directions. The last line in Eq. (|7l[) is obtained by an application of Wick's theorem 
|64| to the Gaussian expectation values {a^wf)w The 'propagators' IL, play the role of 



generalized diffusion poles. They are defined through 

-Dd 2 x + 2ie(2cos#(0) -3cos0(x))] U^x^) = 8(x - y), 

-Ddl - 2ie cos 6(x)] U 2 {x, y) = 8(x - y). (72) 

Without going into details we remark that the relative minus sign between IT and II2 in (|71~D 
derives from the different symmetries of the matrices Wi and w 2 under matrix transposition 
(cf. Eq. ©). 

We will not proceed any further analytically. In order to quantitatively evaluate the 
Il-dependent corrections to the DoS, one would have to compute the generalized diffusion 
poles (|72|). Due to the presence of the spatially varying terms ~ cos^, a general solution of 
the differential equations is difficult^. Nonetheless, quite a few characteristic properties of 
the DoS corrections can be deduced from ([H]) simply by inspection: 

For asymptotically large energies e ^> E c , one expects no influence of the superconductor 
on the normal metal. Indeed, in that limit, cos 6 — > 1 implying that a) the Usadel DoS 
becomes metallic and b) III — II2 — > 0, i.e. no quantum corrections to the DoS. 

For intermediate energies just above the minigap edge E g , cos 6 varies smoothly as a 
function of position. In this regime Eq. ( [7TD gives corrections of 0(g~ l ) to the DoS whose 
quantitative evaluation is difficult. 

Finally, let us consider subgap energies, e < E g (remaining of course outside the regime 
e <C E c ). Here, according to the quasiclassical analysis, the DoS vanishes, implying that 
cos# is purely imaginary and the effective action ( [7DD purely real. As a consequence, the 
propagators IT are real, too, and the DoS, as computed according to ([H]), vanishes iden- 
tically below the quasiclassical minigap edge. In other words, the perturbative inclusion 
of first order quantum corrections does not give rise to the appearance of states below the 
quasiclassical edge. The vanishing of Re (cos 6) actually suffices to demonstrate that the 
robustness of the gap pertains to all orders of perturbation theory. 

This conclusion presents something of a puzzle: Taking into account quantum correc- 
tions, the above result indicates that the minigap remains fixed at energy E g . The latter is 
determined by the bare value of the diffusion constant, D = vpljd. However, the intuition 
afforded by the one-parameter scaling theory of localization |73[ suggests that observables 
such as the tunneling DoS should depend only on the value of D renormalized by weak 
localization corrections. In bulk normal metallic samples, weak localization corrections (to 



13 In fact, relatively standard techniques [75| may be employed for the solution of Eq. (^). This 
follows from the fact that, upon substitution of the Usadel solution Eq. ( p!2| ), these equations 
are classifed as 'Lame' equations [76]. In comparatively simple situations, such as asymptotically 
large energies e 3> E c , infinite SN- rather than finite SN-systems, and so on, analytical solutions 
are available. However, in order not to diversify the discussion unnecessarily we do not elaborate 
on these cases. 
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two-particle properties) stem from mechanisms of quantum interference between trajecto- 
ries connected through a time reversal operation, (see fig. |i"6|(fl)). In the present case, weak 
localization effects can arise due to the interference of particles and holes (see fig. [T6|( b)). 



Since holes bear similarity with time reversed particles, there is no conceptual difference 
to the above N-interference mechanism, and one expects a standard renormalization of the 
diffusion constant (albeit already on the level of one-particle properties). Yet, according to 
the analysis above, the minigap edge, a function of the unrenormalized D, is robust against 
perturbative quantum corrections in the particle/hole channel. 





FIG. 16. Renormalisation of the diffusion constant, D, by interference of (a) trajectories with 
their time-reversed counterparts, and (b) particles and holes. 

In fact, the absence of weak localization corrections to the minigap edge signals the failure 
of the perturbation theory. To properly identify quantum weak localization corrections to 
the diffusion constant, and therefore the minigap edge, it is necessary to renormalize the 
saddle-point equation itself. This situation parallels that encountered in the study of the 
renormalization of the gap in a dirty bulk superconductor where quantum corrections (in 
the Cooper channel) lead to a renormalization of the gap equation (see e.g. Ref. f77j). In 
the present case, weak localization corrections to the minigap edge are obtained within a 
renormalization group procedure. Since, operationally, this procedure is somewhat technical, 
its description has been made the subject of appendix [E|. 

The renormalization group procedure described in appendix [E] may be employed safely 
down to energy scales in excess of E c . However, at energy scales in the vicinity of the 
minigap, the Cooperon propagator depends sensitively on the geometry and it becomes 
necessary to include the additional flow in R (the matrix rotating to the RG-affected saddle 
point of the theory), coupled to that in D, as the cutoff is lowered towards E g . Although the 
manner in which such renormalisation processes are included self-consistently lies beyond 
the scope of this paper, the outcome of the RG procedure can be summarized as follows: 
Treating quantum interference correction within a RG-scheme leads to a shift of the minigap 
edge. The overall structure of the gap edge (e.g. the non-analytic behaviour of the DoS 
in the vicinity of the DoS) is maintained. In particular, no states are found below the 
(renormalized) gap edge. 

These findings leave us with the question whether indeed, no states exist below the 
(renormalized) gap edge or whether the computation simply has not been accurate enough. 
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Although a quantitative analysis of this question is beyond the scope of this paper, we 
believe that the second option is the correct one: To find states below the minigap, one 
must account for contributions to the action which cannot be accessed by a perturbative 
shift of the inhomogeneous saddle-point. Contributions of this kind have been identified 
in bulk normal conductors as soliton-like configurations of the Q-matrix fields, and have 
been associated with a rare class of states which are described as "anomalously" or "nearly" 



localized within the metallic phase 



Poorly contacted to the superconductor, 



these states are able to exist at energies below the minigap and generate contributions to 
the average DoS exponentially small in g. Although we see the SN system as a useful and 
challenging arena in which to investigate the localization properties of such rare states, their 
consideration lies beyond the scope of this paper. 




prelocalized states 



C-mode dominant 



exotic' diffusion modes 



weak localization type corrections to DoS 
> 



E 9* 



FIG. 17. Schematic indication of the role of the various DoS corrections. E g * denotes the 
renormalized gap edge. 



Before leaving this section, let us in summary list the - admittedly diverse - set of 
a-type fluctuation mechanisms renormalizing the single-particle properties of mesoscopic 
SN-structures (see figure [171) : 



For energies e ^> E c , the N-component behaviour is largely metallic. However, the 
presence of the superconductor is exerted in terms of massive quantum corrections 
to the DoS and other single particle properties. The larger e is, the smaller are the 
corrections. 

Energies just above the minigap edge E g are the most difficult to analyse. Quantum 
corrections to quasiclassics are carried by diffusion type modes which - due to the 
pronounced energetic and spatial inhomogeneity of the DoS - are difficult to treat 
analytically. By perturbatively including such corrections, one obtains corrections to 
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the DoS above the gap. Both the position of the gap and the vanishing of the subgap 
density of states remain unchanged. 



• By embedding the diffusive modes into an iteration of RG-analyses and solutions of 
renormalized Usadel-type mean field equations, one arrives at a shifted minigap edge. 
The non-analyticity of the gap is maintained - that is, there are no smooth DoS 'tails' 
leaking downwards out of the sharp edge. 

• Presumably, 'nearly localized' subgap states can be found with a probability that is 
exponentially small in the metallic conductance g. 

• Eventually, for energies e <C E c , the fluctuation physics is governed by the C-mode 
whose impact on various physical observables (for non-proximity effect SN-structures) 
has already been discussed in the literature. 

It is important to question whether the above corrections can be made experimentally visible. 
As far as the DoS is concerned, the answer must be a conservative one: the chances are that it 
will be impossible to separate the high energy 1/o-corrections from the Usadel background. 
Furthermore, for good metals (g ^> 1), finding nearly localized subgap states will also 



be difficult, since, as shown in Refs. pOJS^JSS], disorder configurations leading to nearly 



localized states are exponentially rare. Thus, as far as the mean DoS is concerned, the 
above fluctuation contributions will probably be hard to detect. However, the primary 
purpose of this section has been to demonstrate that a variety of interference mechanisms 
adding to the standard quasiclassical picture exist in principle. If and to what extent these 
fluctuations give rise to observable changes in single particle properties other than the DoS 
(e.g. the Josephson coupling characteristics) represents a subject of future research. 

We now leave the issue of the renormalization of single particle properties and turn to 
the discussion of correlations between more than one Green function, as described by the b- 
and c-type fluctuations. 



B. b-Type Fluctuations: The Goldstone Mode 

In this section we discuss the class of fluctuations around the Usadel saddle point which 
above has been denoted by 'type b)'. Unlike the a-fluctuations, fluctuations of type b) induce 
correlations between different Green functions. What makes the b-fluctuations particularly 
important is their Goldstone mode character: In the limit of vanishing energy difference 
between the considered Green functions, these modes become truly massless, a signature for 
the presence of pronounced mesoscopic fluctuations. 

Consider the action fl33f ) in the simple case u + = A = 0. Obviously, any transformation 
Q — > TqQTq 1 leaves the action invariant provided that [T, of h ] = 0. Among the group of 
matrices G = {T = T <g> II ph }, there is a subgroup H C G , [H , o~f] = which not only 
leaves the action invariant, but also the saddle point solution Q itself. As a consequence, 
fluctuation matrices contained in the subgroup H are completely ineffective and do not 
couple to the theory. However, the elements T of the coset space Go/H do generate non- 
trivial transformations of the diagonal saddle point. Moreover, in the limit To = const., 
these transformations do not alter the action - they are Goldstone modes. 
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Being Goldstone modes, the effective action of the To's can only contain gradient terms 
and mass terms induced by sources of symmetry breaking, such as finite u + and A. The 
actual structure of the action depends crucially on its behaviour under time reversal. For 
the sake of simplicity, we focus here on the two pure symmetry cases: 

i) Orthogonal symmetry: The action is time reversal invariant, A = 0. 

ii) Unitary symmetry: Time reversal invariance is broken, |A|L$q 1 ^> g^ 1 ^ 2 . 

Here, $o is the flux quantum. Note that A denotes the vector potential with account 
for the phase difference between the superconducting terminals, so that phase differences 
A(/)$q 1 ^> g~ x l 2 suffice to drive the system into the unitary symmetry class. The reason that 
g -1 / 2 appears as a measure for the strength of the perturbation is that, for | A|L$q 1 3> g^ 1 ^ 2 , 
the dimensionless coupling constant of the symmetry breaking operator in the action exceeds 
unity (see Appendix |B|). Alternatively, one may say that under these conditions, the mass 
of the 'Cooperon' greatly exceeds the level spacing. 

The derivation of the effective action, So, of the Goldstone modes is somewhat technical 
and has been deferred to Appendix [D|. Here we merely state the result, 

7T f r "I 

S [Qo) = -c S ym- J [v n D str (dQ dQ ) + 2iu + u str (Q ag)\ , (73) 

where D = ^(1 — q + - q— ) plays the role of a space dependent diffusion coefficient. The 
variation of D with both position and energy is shown in fig. |Tj| for an SN junction with a 
typical choice of material parameter. In addition, v is the space dependent, local DoS, as 
displayed already in fig. ||. (Notice that in the bulk of S both the diffusion constant and v 
vanish. Hence, the support of the action of Qo is restricted to the N region.) 




FIG. 18. The effective diffusion constant, D, for the Goldstone modes in an SN junction as a 
function of both energy and position, for 7 = 0.1. 
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Further, in the case of 



i) Orthogonal symmetry: c sym = 1, the matrices Tq are eight dimensional and obey the 



time reversal symmetry relation (50), and of 



ii) Unitary symmetry: c sym = 2, the matrices T are four dimensional (i.e. they do not 
carry a tr-index structure) and Eq. ( |50"D is meaningless. 

In either case the matrices T obey the restricted version of the symmetry relation (pS|)p^| 

Tl = Vo T Q %\ (74) 

where 



770 = E\[ ® <x 3 ar + E. 



22' 



In summary, we see that the symmetry of the Goldstone fields is identical to those of the 



standard Efetov Q-matrix manifolds [63j. In other words, by freezing out the ph-degrees of 
freedom, the large 16-dimensional cx-model manifold collapses to smaller ones of dimension- 
ality 8 (4) which are symmetrically identical to those encountered in orthogonal (unitary) 
applications of the standard cr-model. 

Besides the general symmetry relations, a further condition to be imposed on the fields 
is that they obey Neumann boundary conditions d±T (x) = at all SN- interfaces. The 
derivation of these boundary conditions is discussed in Appendix [B]. 

So far the discussion has been for a general SN geometry. In order to actually demonstrate 
how the b-modes generate mesoscopic fluctuations we next consider a specific example, 
namely, the problem of DoS fluctuations above the minigap edge in an SNS-structure. 



1. Level-Statistics in SNS- structures 



It is well known that the single particle spectrum of mesoscopic, purely normal systems 
is governed by various types mesoscopic fluctuations (see e.g. Ref. |3| for review). The 
fluctuation behaviour can be characterized conveniently in terms of correlations ~ (u(e — 
uj/2)v{e + ui/2)) between the DoS's at different energies. Extensive analyses of correlation 
functions of this type have shown that the DoS-correlations become increasingly universal 
in character, the lower the energy separation u (a fact that follows heuristically from the 
interpretation of to as an inverse time scale). In particular, for energies e < E c the correlations 
become fully universal in the sense that they depend on nothing more than the mean (and 



14 At first sight it seems like we are facing a problem here: A rotation matrix cannot be of the 
form T = Tq (g) IL ph and simultaneously obey the general relation d49|). The reason that matrices of 
b) type are nonetheless permitted is that the relation ( [49|) is in fact too strict. What matters is the 
restriction of the symmetry to the boson-boson and the fermion-fermion block of the matrices T 
(see the corresponding discussion in Ref. |fT3[ .) With regard to these bf-diagonal blocks, matrices 
Tq obeying Eq. (74) are compatible with Eq. 
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constant) level spacing d and the fundamental symmetries of the system. This is the regime 
of Wigner-Dyson statistics. For larger energies, the Wigner-Dyson behaviour crosses over to 
other and less universal types of statistics. Nevertheless, the correlations remain energetically 
long-ranged in the sense that they decay algebraically as a function of u>. 

Here we ask to which extend this behaviour carries over to the fluctuation behaviour of 
the SNS DoS in the vicinity of the minigap. As compared to normal metals, the situation is 
more intricate in that already the mean DoS is affected by mechanisms of quantum coher- 
ence. A conceivable situation is that the (strong) modes of quantum interference, giving rise 
to the particular structure in the mean DoS, decouple entirely from the modes responsible 
for DoS fluctuations. Another possibility is that one might end up with some kind of in- 
separable conglomerate of modes of interference, and thereby fundamentally non-universal 
types of spectral statistics. Here we demonstrate that the true picture lies somewhere 
halfway between these two extreme options: It is still possible to identify a regime of uni- 
versal Wigner-Dyson statistics, albeit superimposed on an energetically non- uniform DoS. 
However, its range of validity shrinks down to a small energy window beyond which the 
correlations do become entirely non- universal. 

To obtain specific information about spectral fluctuations in an SNS-structure, we apply 
the action ([T3|) to the analysis of the SNS geometry, shown in fig. £| The set of field config- 
urations obeying Neumann boundary conditions at the SN-interfaces obviously contains a 
subset with trivial spatial dependence, T = const. - the 'zero mode'. The zero mode action 
reads: 

S [Q = const.] = -i||!±_stro[Qo<]- ( 75 ) 

Here d(e) = (J z/(e)) _1 denotes the average level spacing, which, in contrast to the purely 
normal case, is now energy dependent. Note the similarity with the action, (EO), of the 
C-zero mode. The difference is in the pre-factor and in the physical spaces in which the 
matrices Q and Q° C) respectively, operate. 

In order to demonstrate any significance of the isolated zero mode action, we have to show 
that it is separated by an energy gap from the action of all other field configurations. At this 
point, the boundary conditions begin to play a crucial role. Expanding the fields in terms of 
cosines, that is, a complete set of functions compatible with Neumann boundary conditions, 
we see that, next to the zero mode, the field with least curvature varies as ~ cos(27rx/L). 
Due to the presence of the gradient term, the action associated with this configuration is of 
0(D/ \L 2 d). Thus, for u + <C E c , the zero mode action is energetically gapped against all 
fluctuating contributions and plays a dominant role. 

We will see in the next section that, due to the presence of the c-type fluctuations, the 
range of stability of the zero mode is actually much smaller than uo + < E c . Yet, restricting 
ourselves for a moment to the consideration of the zero mode action, we can, without any 
further calculation, draw immediate conclusions about the level statistics in SNS-systems 
over small correlation intervals, u + . In fact, actions of the structure ( [75] ) are standard in 
applications of the cr-model in N-mesoscopic physics: they appear a) whenever a model may 
be subjected to a zero mode approximation or b) when one is dealing with cx-model analyses 
of a single random matrix ensemble. With regard to spectral statistics, the existence of the 
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zero mode action (]75|) implies that level correlations for small energies are of Wigner- Dyson 
type. 

Furthermore, a comparison of the action ( fT5| ) with the analogous action for N-systems 
|64| shows the correlations to depend on an average level spacing that is effectively halved. 
This is due to the strong 'hybridisation' of levels at energies ~ ep ± e induced by Andreev 
scattering at the SN-interface. 

A more comprehensive discussion of level fluctuations, including the differences to the 
types of spectral statistics found in N-materials, will be given after the c-type fluctuations 
have been incorporated in our analysis. 



C. c-Type Fluctuations: Quantum Corrections to Level Statistics 



In the previous section, an effective action for b-type fluctuations was derived, and the 
latter were shown to be Goldstone modes of the theory. Furthermore, for the SNS geometry 
at energy scales u + <C E c , the effective action was shown to be dominated by a zero mode 
which established universality of level statistics within the ergodic regime. Higher modes 
give rise to non-universal corrections on energy scales u + ~ E c . At the same time, c-type 
fluctuations, that is, fluctuations that commute with neither cr ph nor cr ar ), also become 
important. 

The aim of this section is to examine the role played by c-type fluctuations in limiting 
the regime over which universal correlations persist^. Since, for states above the minigap, 
c-type fluctuations incur a mass which is of order e/d > g ^> 1, fluctuation corrections to the 
universal level statistics can be treated within a perturbative manner. Our approach will 
be based on the perturbative treatment introduced by Kravtsov and Mirlin |79| in studying 
similar corrections in normal disordered conductors. 

We begin by employing the general parameterization 



Q = T Tq • a 



phrp—lrp — l 



T 



cxp 



fj,=0 



(76) 



where Qo = Toq ■ a^Tg" 1 represents the spatially homogeneous zero mode. Here we have 
applied the notation o-q = 11, and, separating the zero mode, we impose the requirement 
that the fluctuations obey the constraint / W = 0. Here, as before, q = q ■ cr ph represents 
the saddle-point, or Usadel, solution. 

Applying this parameterization, and expanding to quadratic order in the fields W, the 
total correction to the zero mode action (1751) takes the form 



It is conceivable that the unfolding procedure of the previous section, that allows for a energet- 
ically inhomogeneous mean DoS, may itself be a source of decorrelation on energy scales u;+ ~ E c . 
Such a mechanism for non-universal corrections to the level statistics would be separate to that 
described here for the c-type fluctuations and so is not contained within our present analysis. 
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5S[Q ,W] 



7TU 



str{D {[dW, q] 2 + 2[dW, W]qdq) 
-2itq [W, [W, af] 
+4iuj + U af[W,q}} + O(W 3 ), 



(77) 



where Uq = T^~ 1 a^ r T . In contrast to a normal conductor, inhomogeneity of the saddle-point 
solution q allows a term linear in W to survive in the action. The presence of this term has 
important consequences on the range over which level correlations are universal. 

To proceed, it is convenient to further separate c-type fluctuations into two classes, 
W = W A + W D : 



Modes diagonal in ar-space ([W / / f,o"3 

rD _ urD ^ a Ph + W D 3 



0), but off-diagonal in ph-space, are termed 



D modes; W u = Wf 



ph 



Modes off-diagonal in ar-space ([W A , uf]+ = 0), but diagonal in ph-space, are termed 



A modesQ W J 



W A 



l ph + W A ® at- 



On the level of the quadratic action, no mixing between these modes occurs. Finally, con- 
siderations analogous to those presented in connection with the a-type fluctuations show 
the spectrum of these modes to be discretely spaced, where the typical 'level distance' is 
of 0(E C ). In passing we note that the c-modes, as introduced above, are in fact not com- 
plementary to the a-modes. For example, the above W D, s contain modes oc IL ar , which, 
by definition, belong to type a). However, the present analysis, regarding the impact of 
non-universal fluctuations on the b-type Goldstone mode, does not require a separation of 
the a- and c-modes and it is sufficient to continue with the present definition of the c-modes. 

Once again, to keep our discussion simple, we limit consideration to pure symmetry 
classes of either orthogonal or unitary type. This leads to a simplification of the effective 
action (|77|) allowing an explicit integration over the fluctuations W. Specifically, for pure 
symmetry, the action takes the form 5S = Sa + Sd, where 



Sa 
Sd 



7TZ/ n 



str 
str 



W A (-dbd) W A + tuj + W A a^([q 3 } + - [q 3 }_)U 



W° [—Dd z - 2ieq 3 ) + 2u+W"q 1 U [ 



rD 



rD, 



(78) 
(79) 



and D is the space and energy dependent diffusion constant that has been introduced in 
section |VII B| for the Goldstone modes. Note that the W 3 fluctuations do not couple linearly 
to Uq in the pure symmetry case and so may be dropped. 

As can be seen from the general structure of the action, c-type fluctuations in the vicinity 
of the minigap are generally 'massive', that is, governed by an action which is at least of 



The denotation 'A' respectively 'D' modes is again motivated by Cartan's classification scheme 
of symmetric spaces (see Ref. for a discussion of the scheme focusing on its application to the 
symmetry classification of SN-systems.) 
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order e/d ^ g 3> 1. It is thus permissible to treat these fluctuations in a simple Gaussian 
approximation. Applying the shift operations 



K - < - ^Po(-dDd)-^\[q 3 ] + - [q 3 ]_)U , (80) 

W 2 D -> W 2 D - ^ + P 2 (-^ 2 )- 1 g 1 f/ , (81) 

where P M represents a projector onto the field space of then performing the Gaussian 
integral, we obtain the renormalized zero mode action 




S[Q ] = S [Q ] str[<, Q ] 2 , (82) 



where «(e) ~ 0(1) represents a constant that depends on the sample geometry. 

Eq. (|82"D has a structure equivalent to that found in the study of universal parametric 
correlation functions and explicit expressions for the two-point correlator of DoS fluctuations 
for both orthogonal and unitary ensembles can be deduced from Ref. Qualitatively, 
the additional contribution in (SJ3) counteracts the zero-mode fluctuations for non-vanishing 



frequencies u + . 

Furthermore, we find a marked difference in the manner in which level correlations are 
suppressed as compared to the purely N case. Already for energy separations u)+/d(e) ~ y/g, 
the zero-mode integration is largely suppressed which manifests in an exponential vanishing 
of the level correlations on these scales. This represents a qualitatively smaller energy scale 
than that for the purely N case, for which Wigner-Dyson statistics prevail all the way up to 
frequencies u + ~ E c . In addition, the exponential suppression of correlations in the SN case 
differs from the purely N case, for which the Wigner-Dyson statistics are succeeded by other 
forms of algebraically decaying spectral statistics in the high frequency domain u + > E c 
|85| . Note that a similar phenomenon of zero- mode suppression has recently been observed 



by Skvortsov et al. |52) in their analysis of the level statistics of normal core excitations in 



type II superconductor vortices. 



VIII. DISCUSSION 

In conclusion a general framework has been developed in which the interplay of meso- 
scopic quantum coherence phenomena and the proximity effect can be explored. The con- 
nection between the conventional quasiclassical approach and the field theoretic approach 
adopted here has been emphasized. In applying the effective action we have introduced a 
classification of different modes of fluctuations. 

To keep our analysis simple we have focussed on a regime in which the contact between 
the superconductor and normal regions is metallic, and where A ^> E c . Experimental 



17 With reference to the specific correlation function i?2, we remark that only massive fluctuations 
in the ph-sector contribute to connected correlators of the form {G^Q^j allowing such terms to 
be neglected. 
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analyses are often carried out in the complementary regime where tunnel barriers separate 
S and N and/or A <j E c . It is straightforward to modify the theory so as to accommodate 
tunnel barriers, small order parameters and, in fact, altogether different sample layouts. 
However, in order not to diversify the present exposition of the formalism even further, we 
have restricted ourselves to the analysis of the relatively simple systems discussed above. 
Whereas certain of our conclusions (e.g. the existence of a Wigner-Dyson regime of spectral 
correlations) carry over to the case of barrier separated SN systems, others do not. More 
specifically, reducing the strength of the order parameter below E c affects both the behaviour 
of certain of the fluctuation classes discussed above and the spatial structure of the solutions 



of the Usadel equation |36| . Rather than attempting to set up a most general 'phase diagram' 
of mean field and fluctuation regimes - given the diversity of SN-systems with qualitatively 
different physically behaviour, certainly a fruitless task - it is more efficient to treat different 
problems individually, i.e. to start out from the most general form of the action (^) and to 
restrict the analysis to those fluctuation modes that encompass the physics particular to the 



problem under consideration (see, e.g. Refs. [|89| and for recent examples.) Whether or 



not a certain type of fluctuation around the Usadel saddle point is 'relevant' or not can be 
deduced from the way it couples to the different contributions to the action. 

Finally, we note that, in this paper we have focussed on the influence of mesoscopic 
fluctuations on the proximity effect in disordered SN-structures. However, mechanisms of 
quantum interference analogous to those discussed here also induce mesoscopic fluctua- 
tions in irregular clean or "quantum chaotic" structures. Moreover, the proximity effect is 
strongly influenced by such coherence phenomena allowing them to be employed as a poten- 
tial probe of chaotic behaviour p7|j90|1 . Can the framework developed above for disordered 



SN-structures be generalized to account for chaotic or ballistic SN-structures? To address 
this question we should begin by recalling the properties of normal chaotic structures. 

In fact the connection between the statistical field theory of normal disordered conductors 
and ballistic chaotic structures was motivated by the quasiclassical approach of Eilenberger 
discussed previously. Recognizing that the Usadel equation could be associated with the 
equation of motion corresponding to the saddle-point of the action of the diffusive non- 
linear a-model | 8*0|] , Muzykantskii and Khmel'nitskii proposed that the Eilenberger equation 



could be identified with a ballistic analogue of the non-linear a- model action [p4 j. In this 
case, the diffusive character of the action was replaced by a kinetic operator. Their work 
found support in subsequent investigations based on the study of energy averaged properties 
of (again normal) chaotic structures which led to a microscopic derivation of the ballistic 



action |87[] . Taken together, these studies showed that, while density relaxation in disordered 
conductors is diffusive, in general chaotic structures it is governed by modes of the irreversible 
classical evolution operator. 

The generalization of the ballistic field theory to encompass the proximity effect follows 
naturally from the ideas presented in this paper. Expanding the field space of Q to ac- 
commodate particle/hole degrees of freedom and, as in the diffusive model, introducing the 
inhomogeneous order parameter A, the effective ballistic action takes the form 



S[Q] =i^j str [ 2iT < ® °fi H , T- 1 } - of ® (A + e + Q 



(83) 



where {H, } represents the Poisson bracket of the classical Hamiltonian, and the supermatrix 
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fields Q = Taf ® af^T" 1 depend on the 2d — 1 phase space coordinates parameterizing the 
constant energy shell, xn = (r,p) 2 d-i- (With this definition the integration measure is 
normalized such that J = J dx» — 1.) In the presence of a Gaussian distributed 5-correlated 
impurity potential, the ballistic action is supplemented by a further term corresponding to 
a collision integral [[M 



Sou = ^ J ^~sT str [g(n ' r)g(n ' r)] ' (84) 



where n = p/|p| and Sd — f dn. Indeed, for a strong enough impurity potential, dr <^ 1, 
a moment expansion of the action recovers the diffusive action. Varying the action with 
respect to Q, and applying the identification g(n, r) <-» Q(n, r), the saddle-point equation 
of motion coincides with the Eilenberger equation of transport, Eq. @. 

Although, in principle the ballistic action represents a complete theory of statistical cor- 
relations in chaotic SN-structures, an analytical description of the modes of the classical 
evolution operator has proved difficult to construct. In particular, the sensitivity of weak 
localization corrections to mechanisms of "quantum diffraction" and "irreversibility" in nor- 
mal clean chaotic structures has proved difficult to quantify f40|| . In the SN-geometry the 



same mechanisms have a dramatic effect on the single-particle properties of the device such 



as the minigap structure in the local DoS |37] ( see the discussion in section [II B|) . For this 
reason, we believe that SN-structures may provide a versatile arena in which properties 
quantum chaotic systems can be explored. 
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APPENDIX A: BOUNDARY CONDITIONS OF THE USADEL EQUATION 

In addition to the transport equations provided by quasiclassics, it is necessary to specify 
boundary conditions at the SN-interfaces. For the Eilenberger equation, at a planar SN- 
interface with an arbitrary transmission coefficient, T, these conditions have been derived 
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by Zaitsev |^7|. Note that these boundary conditions cannot be obtained using the standard 



quasiclassical Green function (y) alone. Instead, one must go back to a more microscopic 
formulation. 

Following the general philosophy of this section we shall not review the (somewhat tech- 
nical) derivation of the boundary conditions but merely formulate the main results. The 
Eilenberger Green function, g(n, r), may be separated into symmetric and antisymmetric 
parts, g = g s + g a , with respect to the operation vp — ► -vp. The antisymmetric part, g a , is 
continuous across the interface. In passing we note that this results in the conservation of 
the supercurrent density, 

j = -f|(ntr ( T 3 P V(n,r)) n , (Al) 

across the interface. In contrast, the symmetric part experiences a jump depending on the 
transmission coefficient. The resulting conditions are: 

9a(+) = g a (-) = g a , (A2a) 
g a - g a g a ) + j(g,(+) - <? S R) 2 } = j [&(-), g.{+)], (A2b) 

where R = 1 — T is the reflection coefficient, the r,a superscripts have been dropped, and 
<7a(±) denotes the Green function infinitesimally to the left respectively right of the junction. 

For a perfectly transparent (T = 1) interface, both parts of g are continuous. In the low 
transparency limit, T C 1, we have g a ~ T and ( |A2b|) reduces to 

9a = ^[9s(-),9s(+)]. (A3) 

The above boundary conditions simplify further in the dirty limit. As shown by Kuprianov 
and Lukichev [[51]], the reformulation of ( |A2a|) and ( |A3| ) in terms of the Usadel Green function 
leads to the pair of conditions, 

o-(-)gad r g (-) = o-(+)g d r g (+) (A4a) 

T = 1 ^ 1 bto(+),flb(-)], (A4b) 

where cx(±) is the metallic conductance on either side of the interface, and 

e 2 v n v F f 1 T 

Gt = / a(cos a)— cos a, (A5J 

2 Jo R 

is the tunnel conductance of the junction, a is the angle between n and r, and v n is the 
bulk, normal metallic DoS. Note that T may depend here on a. Note also that the second 
condition of Eq. ( |10b|) applies only in the limit T < 1. Lambert et al. 0] have recently 
examined this restriction and how it may be relaxed. In the following, we exclusively consider 
the opposite case of perfect transmittance, T = 1, for which the second condition, Eq. (|10b|) . 
should be replaced by 

9o(+)=9o(-)- (A6) 
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APPENDIX B: SADDLE POINTS AND ANALYTIC CONTINUATION 



This appendix is devoted to a discussion of the question of how the Q-saddle-point con- 
figurations that appear in applications with superconductivity may be accessed from the 
starting-point of the diagonal saddle-point configuration, af <8> cr| r , characteristic for bulk 
metallic phases. In order to specify what we mean by 'access', we first have to summa- 
rize some facts regarding the structure of the field manifold of the u-model. In the polar 



representation of Efetov pl |, a general Q-matrix is parameterized as 

Q = T(a 3 ph ®<)T- 1 , (Bl) 

where af <g> a™ is the generalization of the standard matrix, A = af, to applications with 
a ph-substructure. The rotation matrices obey T G G/H, where G is a group of matrices 
that fulfill various symmetry conditions, dictated by both the internal symmetries of the 
model under consideration and convergence criteria. The group H C G is determined by 
the condition [H, of 1 ® = 0. The most important constraint for the present discussion 
has to do with convergence and reads, 

T f = r{T- x rf x , (B2) 

where 

V = E\{ ® af ® of + E b 2 l 

The functional integration / dQ extends over the coset space G/H. A key question to 
address is whether or not all stationary phase points of the action are accessible within the 
integration domain specified by (|B1|). 

To analyse this issue, we first focus on the simple case of a bulk superconductor in the 
regime |A| 3> e. In this case, the saddle point is unique and reads (cf. section |V B| ) 

Q = af. 

The above question reduces in this case to whether or not there exists a solution to the 
equation, 

af g, 1 ar = T(af ® Or 1 (B3) 

Strictly speaking, no such solution exists. To demonstrate this point, we first note that, since 
the rhs of Eq. ( |B3| ) is trivial in ar, bf and tr-space, it is sufficient to focus on each sector of 
these spaces separately. Without loss of generality, we focus on the retarded-retarded (rr) 
sector, where Eq. ( P3|) takes the form, 

af = TafT- 1 (B4) 

and we have, for reasons of notational simplicity, denoted the rr-restricted rotation matrices 
again by T. The analogous equation for the advanced-advanced block carries an overall 
minus sign on the rhs. Specializing the discussion further to the the fermion-fermion (ff) 
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sector, we find that, in this sector, no problem arises in finding a solution to Eq. ([BJ) with 
the correct symmetries. When restricted to the ff-block, the symmetry relation (|B2| ) takes 
the form 

T't 

A S — 1 S i 

i.e. the ff-matrices are unitary. At the same time, the restriction of Eq. flB~3| ) to the ff-sector 
reads 

which may be solved by a unitary rotation matrix. For future reference, we explicitly write 
the solution as 

T ff = exp(-t6 s af) n e G/H. 

"ff=4 

In the bb-sector, the situation is more problematic: we fail to find a solution to ( [BID with the 
correct symmetries. The restriction of the symmetry criterion ( |B2| ) on the retarded/retarded 
block reads 

1 bb — °3 J bb °3 ) 

which fails to include any T bb that further fulfill 

°1 ~~ J bb°3 J bb • 

The resolution of this problem is provided by analytic continuation. In the derivation of 
the cr-model, the symmetry condition ( P2|) is enforced by convergence requirements. As 
long as no singularities are encountered, the condition may be relaxed, in the sense that the 
integration contours may be analytically continued to regions where the symmetry criterion 
is no longer fulfilled. Supposing now we are integrating over the subset (cf. fig. [19]) 



{T = expH# bb( x 2 ph )|# bb G ill} C G/H, 

the saddle point we wish to access is reached by distorting the integration contour so as to 
cross the point # bb = tt/4. In lifting the integration path off the imaginary axis, no singular- 
ities are encountered. Moreover, a closer analysis shows that, in accordance with the basic 
conditions to be imposed on saddle point integrals, the direction of steepest descent is par- 
allel to the imaginary axis. Analogous arguments may be applied to the advanced/advanced 
sector. Consequently we conclude that q does represent a proper saddle point of the Q- 
integration. 
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deformed bosonic contour 
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FIG. 19. Visualization of the deformed integration contour. 



The question remains as to what happens if we encounter more complicated saddle- 
point configurations, such as those with a finite value of e/A or even with spatial variation. 
Although no mathematical proof has been given, we expect that by analytic continuation 
such configurations remain accessible. Independent evidence for the validity of such an 
assumption is provided by physical criteria: The disordered mean DoS may be calculated 
within the framework of the quasiclassical approach, while the solution to the quasiclassical 
equations coincides with the solution of the a-model mean field equations. Suppose now 
that the solution were inaccessible in either block, bb or ff, or both. In this case the 
functional integral would exhibit supersymmetry-breaking on the mean field level and it 
would be obscure how to reproduce the correct results of quasiclassics. Given such sources 
of evidence, we adopt a pragmatic point of view and take for granted the accessibility of the 
Usadel-saddle points. 



In this section we provide explicit solutions of the Usadel equation for two simple, quasi- 
1D geometries: an SN junction and an SNS-junction with coincident phases of the order 
parameters of the S regions (Ay? = 0). 



We begin with the SN geometry. As with the case of the bulk superconductor, we 
introduce an angular parameterization for q(x): 



APPENDIX C: SOLUTIONS TO THE USADEL EQUATION 



1. SN junction 
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q{x) 



(sm9(x), 0, cos9(x)). 



(CI) 



Since in this geometry the phase of the order parameter is spatially constant, the gauge 
transformation of section [VI] allows us to set the second component of q to zero. The 
saddle-point equation for q, Eq. ([HP , becomes a sine-Gordon equation, 



(C2) 



^-d 2 x 9 + ie sin 9 + A(x) cos 9 = 0, 

while the boundary condition at the interface, Eq. ( |45b|) , becomes 

a s d x 9(0-)=a n d x 9(0 + ). 
Also, the symmetry relation, Eq. (^), becomes 

9Jx) =7T- (0+Y. 



(C3) 



(C4) 



In addition, there are further boundary conditions at infinity, at which the bulk values of 
the angle are approached, so that 



9{x) 



9 S , x — 

9ni 



-oo, 
oo, 



(C5) 



where 9 S is defined by the equations for the bulk order parameter, eqns. (p7|), and (#„,)+- = 
0, 7r. Note that, although the conditions ( |C5|) takes different forms in the ar-sectors, the 
existence of the relation ( |C4| ) means that we need solve only for the retarded component, 
9 + , and in the following we drop the + subscript. 

The solution of the sine-Gordon equation, Eq. (p3), with boundary condition flC5p, is of 
the following (solitonic) form: 



9(x) 



9 S + 4tan" 1 



exp 



D s 



x tan 



9(0)- 



4 tan" 1 



exp 



-2ie, 



x tan 



0(0) 



x < 
x > 



(C6) 



where R = |A| 2 — e 2 as before. 

The integration constant #(0) is fixed by imposing the condition ( p3|) at the interface, 
to give 



sin = 7, 



. 9(0) 
sin — — 



where 7 is a parameter representing the mismatch between the two materials: 

7 



D n /D a 



(C7) 



(C8) 



In the limit 7 — > 0, the bulk value of the angle, 9 S , is imposed asymptotically at the interface, 
9(0) —>■ 9 S - the boundary condition becomes 'rigid'. 

By Eq. (ESf), the local DoS is obtained from the relation, 



u(x) = u n Recos9(x). 



(C9) 
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2. SNS Junction with Coincident Phases 



We turn now to the geometry of an SNS-junction, of width L, and with coincident phases 
of the order parameters in the S regions. Overall symmetry about the origin leads to the 
condition, 

0*0(0) = 0. (CIO) 

Since an identical condition holds at a normal-insulator interface, the solution here also 
applies to an SNI junction of width L/2. There are further conditions at infinity, 

9(x) -> 9 S , \x\ -> oo. (Cll) 

The solution of Eq. flU2|) for 6(x), incorporating eqns. QU10| ) and ( |C11[ ), is as follows: 



s + 4tan- 1 fexp (-*[*&(\x\ - L/2)) tan e W 2 )~ A , |a;| > L/2, 
9{x) = { V V V s / ) (C12) 

2 sin" 1 (sin 9 -f sn (i ^ x + K (sin e -f) , sin 9 -f) ) , |x| < L/2. 

Here and sn are the complete elliptic integral of the first kind and the Jacobi elliptic 



function, respectively (see [78]). The two integration constants, 9(0) and 9(L/2), are related 
by 

and the conditions at the interfaces, corresponding to Eq. flC3|) , give the further relation 

. (9(L/2)-e.) rz-( . 2 e(L/2) . 2 »(o)\ 1/2 

sm , =7 /_/ sm ^ sm2 _l , (C14 ) 



where the parameter 7 is defined by Eq. ( |C8|) as before. The integration constants may then 
be determined by numerical solution of eqns. ( |C13|) and ( |C14|) . 



APPENDIX D: EFFECTIVE ACTION OF THE GOLDSTONE MODE 

The subject of this Appendix is a derivation of the effective action for the Goldstone 
modes, represented by rotation matrices T such that [T , cr ph ] = 0. We consider separately 
the cases of orthogonal and unitary symmetry. 

1. Time Reversal Invariant Action 

We begin with the case of orthogonal symmetry. Substituting the ansatz Q = TqQTq 1 
into the effective action, Eq. (|33D, we notice that we obtain two Tq dependent terms: a 
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gradient term, and a term proportional to uo + . Note that the two remaining vertices, pro- 
portional to e and A, commute through and do not couple to T . Focusing on the gradient 
term first, we obtain 

-nDv n r nDvr, 



-J st T(dQdQ) = -^JstT((d+[T 1 (dT ) J .])Q(d+[T 1 (dT ),.])Q) 

( 2 \ 

([T- l (dT ),Q]) 2 +2 [T 1 (dT ),Q]dQ^ 



txDv, 



str 



+ T —independent . (D 1 ) 



V a) b) J 



We first observe that the contribution b) vanishes, through the use of the condition Q 2 = 11. 
Turning to the a)-term, we write 

a) = str {[T ~\dT ),Q] 2 ) = str (dQdQ) , 

where d is defined to be a derivative that acts only on T . Making use of the ph- 
commutativity of the To's, we next introduce 

Q = ToQT, 1 = ^Aq • a ph Q , 



where 



Qo = Tqg™T q 1 



and Aq = q + — q . Tracing out the ph-indices, we now obtain 

a) = str (dQdQ) = Aq str (dQ dQ ) . 

Substituting this expression into Eq. ([Dl|) and using the fact that Aq • Aq = 2(1 — q + ■ q ^ 
we finally arrive at 

TlDVr, f ,„„„„ s 7TZ/ 



/7TV f ~ 
str (dQdQ) = I Dstr (dQ dQ ) (D2) 

as the gradient term of the Goldstone action. Here, 

D=|(l-q + -q_) (D3) 

plays the role of an effective (and generally space dependent) diffusion coefficient. Note that 
the symmetry relation (pE8| ) implies that this diffusion coefficient is real. 
We next turn to the computation of the co> + -dependent vertex: 



J str (Qaf ® of) = - l I^± J Aq3 stlQ (g 0< ) 

J z/str (Qo<), (D4) 
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where Aq^ is the 3-component of the vector Aq and v = z/„Ag3, the space dependent local 

and (F~ 



DoS (cf. Eq. (|28f)). Combining eqns. 
S [Qo] = - 



we finally obtain 



^Dstr (dQ dQ ) + ^^str (Q < 



(D5) 



as the final expression for the time reversal invariant Goldstone mode action, as included 
in Eq. (j73|) . We emphasize that both the diffusion coefficient, D, and the local DoS, u, are 
space dependent. 



2. Broken Time Reversal Invariance 

If the invariance under time reversal is broken by an external magnetic field and/or 
significant phase differences between the adjacent superconductors, rotations operating non- 
trivially in tr-space are frozen out. More precisely, if the total flux threading the system, 
$/$(), exceeds l/g 1//2 , symmetry breaking contributions to the action with coupling constants 
greater than unity appeaiQ. In the infrared limit, fluctuations coupling to these vertices 
become inessential and only T-matrices with tr-block diagonal structure survive, through 
their commutativity with the tr-symmetry breaking operators: 

T = T x ® El[ +T 2 ® E*. 

The symmetry relation flo*0| ) implies that the two blocks are related to each other by 

Tf = T 2 . (D6) 



Substituting the block diagonal form into the action of Eq. (p3|) , all T- invariance breaking 
operators drop out. Thus, similarly to the preceding subsection, we again arrive at Eq. ( p5| ) 
as an effective Goldstone action. Exploiting the tr-block structure of the T-matrices, the 
action may be written as the sum of two contributions, S = S1+S2, where the subscript refers 
to the tr-index. Due to Eq. (|D6| ) and the invariance of the 'str' under matrix transposition, 
the two contributions are identical and we obtain 



SolQi 



7TU ~ 

— -Dstr (dQodQo) — i7ru; + z/str (QqO's 



(D7) 



as the final result for the T-non-invariant Goldstone mode action, as also included in Eq. 
In Eq. (|JJ7|) , we have used 



18 The fact that is the relevant scale follows simply from gauge invariance: In the presence 

of fields, the gradient operator appearing in (|D2| ) generalizes to the the gauge invariant form 
Dd 2 — > —D{id — eA/c) 2 . For a static field of strength $/Area, the vector potential A ~ <I>/L 
and the 'diamagnetic term' is of order D(eA/c) 2 ~ E c (&/Qq) 2 . Since the basic energy-unit of the 
action is the mean level spacing, d, the dimensionless coupling strength is of 0((&/&o) 2 E c /d) = 
O((<fr/&o) 2 g). For $/$o > 5 -1 ^ 2 the coupling strength exceeds unity. 
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Qo — t qt 1 

where - for the sake of a homogeneous notation - we have denoted T\ again by T and the 
tr-block Qii by Q. Note, however, that the matrix dimension of the fields in ( |D7| ) is twice as 
small as the one in (|D7| ), as the tr- index structure is missing, and that there is no symmetry 
relation such as Eq. (p0|). 



3. Boundary Conditions 

In order to make the gradient terms appearing in eqns. ( p5| ) and ( p7[ ) well defined, 
boundary conditions at all interfaces between the N-region and external regions need to 
be specified. Whereas the boundary conditions to be imposed at interfaces to insulators 
(<9T = 0) and idealized leads (T = II) have been derived previously f64j] , the interfacial 



behaviour at SN-boundaries has so far not been analyzed. Note that the present analysis 
applies to the boundary condition of the Goldstone mode as opposed to the behaviour of a 
single-particle Green function (the ar-diagonal blocks of the Q-matrices) - The latter has 
been summarized already in appendix 0. 

Some insight into the structure of the boundary conditions may be gained from the fact 
that two types of currents across the SN-boundary may be identified: 

• A potentially non-vanishing electric current. (Two elementary charges flow across the 
interface whenever Andreev scattering takes place.) 

• A quasiparticle current that vanishes (even in the case of a nonzero reflection coefficient 
at the interface). The quasiparticle current democratically counts the flow of electrons 
and holes. Since incoming electrons are either reflected or Andreev converted into 
holes, the normal component of the boundary quasiparticle current vanishes. 

Taking into account the fact that the Goldstone mode does not distinguish between particles 
and holes, one may anticipate that the boundary condition reads d±T = 0, corresponding 
to zero quasiparticle current flow. We confirm this supposition below. 

First it is necessary to decide on the location of the NS boundary. Since superconductive 
behaviour penetrates, in the sense of the proximity effect, into the normal region (and vice 
versa), the position of the boundary is to some extent arbitrary and need not necessarily 
coincide with the material boundary, at which the jump in A appears. However within the 
superconductor both the effective diffusion constant and the DoS vanish on a scale of the 
order of £ = (D/A) 1 ^ 2 , which is much smaller than the diffusion length, L e (for e <C A). 
With a type of coarse-graining in mind for which the details of variation over scales of £ in 
the S region becomes irrelevant, we make a simple choice of the physical SN-interface as the 
effective one for the Goldstone action. 

To derive boundary conditions for the Goldstone mode, we employ the method of bound- 
ary Ward-identities, as previously used in Ref. fl9l| . As is usual with Ward identities, the 



scheme is to subject the Q-degrees of freedom to an infinitesimal gauge transformation and 
to exploit the fact that physical expectation values ought to be invariant, whilst the action 
need not. 
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Specifically, we perform the infinitesimal rotation 

Qo - e~ R Q e R ~Q - [R, Q ] , (D8) 

where 

R(x) = ( R _ +{x) R+ ~ {X) ) ® lt r ,bf, P h, (D9) 
and the matrix structure refers to the ar-indices. A straightforward calculation then shows 



that the effective action ( |D5| ) transforms as 

S [Qo] - S [Q Q ] + 6Sk[Q , R] + 5S s/ n[Qo, R], 

5S„[Q , R] = J [™ nJ Dstr (Rd (Q dQ )) + ^^str (R{Q , af\) 

6S m [Q , R] = ~- [ dSD s/N str [RQ d ± Q ] , (D10) 

Z JS/N 

where / s / N dS" denotes a surface integral over the SN-boundary (induced by an integration 
by parts necessary to shuffle all derivatives from R to Qo) and -Ds/n is the local diffusion 
coefficient at the boundary. In order to use ( piO| ) to construct a boundary condition, we 
consider the functional expectation value 

X = <F(Q (v))>Qo, (DH) 

where F may be an arbitrary function of the matrix Qo(y) at point y 6 N. Whereas both the 
action and F(Qo(y)) need not be invariant under the transformation (|D8|) , the expectation 
value, X, must be. Expanding the expression ( pil| ) to first order in R and omitting the 
matrix arguments in the notation, we obtain 

X -> X + (5F + F(55 N + 5S S /n)> 

and hence 

(5F + F(8S N + 5S S/N ))±0. 

Since the action ^S's/n is singular at the boundary, its contribution to the above expression 
must vanish individually, that is, we have to demand 



F(Q (y)) j dSD s/N str [RQ d±Qo]) = 



/S/N 

for any function F. As R is arbitrary, this can be generally true only if 

(F(Q Q (y)) str [RQ (x)d ± Q (x)]) = 0, \/y G N , x G S/N, (D12) 

where i? may be any matrix of the structure ( |D9| ) . In order to transform Eq. ( pi2| ) into 
a more practical form, by which we mean an effective condition to be imposed on the 
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differential operator governing the action, we subject both the action and the expectation 
value in Eq. ( P12 ) to a perturbative expansion. Introducing 

Qo = e w a-e~ w , W = ^ B ) , (D13) 
and expanding the action to lowest order in B, we obtain 

SolQo]^ J str (bU- 1 b)+..., 

where II is a shorthand for the diffusion type operator governing the quadratic action. Note 
that for the present discussion the detailed structure of II is of no concern. Further, choosing 

F(Q (y)) = str (Q (y)E% ® af) , 

Eq. ( pi2|) takes the form 

(F(Q (y)) str [RQ (x)d ± Q (x)]) - ( str (B(y)af) sti {d ± B{x))) g S 

= Vt/GN.iG S/N, 

where (. . -)b,b stands for a functional average with respect to the above quadratic action. 
Computing the expectation value by means of Wick's theorem, we obtain finally 

d ± Tl(x, y) = Q, Vy G N , x e S/N. 

This is the required boundary condition. It implies that the eigenfunctions of the operator IT 
must be drawn from the set of functions obeying Neumann boundary conditions at the S/N 
interface. In order to install this condition generally, we restrict the functional integration to 
the set of field configurations, B, that obey the same boundary condition, d±B = 0. Since 
T = exp(W), where W is given in Eq. ( |D13| ), an alternative formulation reads d±T = 0. 



APPENDIX E: RENORMALIZATION OF THE MINIGAP EDGE 



In section |VII A 2| the role of a-type fluctuations on the single-particle properties was 
investigated within the framework of a perturbative expansion around the inhomogeneous 
saddle-point solution of the Usadel equation. There it was shown that, in the SNS geometry, 
the minigap induced by the proximity effect is not destroyed by quantum fluctuations. 
However, this calculation failed to account for the shift of the minigap edge resulting from the 
quantum renormalization of the diffusion constant due to mechanisms of weak localization. 
This phenomenon is described below. 

To take into account weak localization corrections, we apply a conventional momentum 
shell renormalization group procedure to the effective action as it is detailed e.g. in pi 



Beginning with the parameterization defined in Eq. (p5|) , we factorise rotations T = T>T< 
into fast T> and slow T< degrees of freedom. Here rotations T> (T<) involve spatial fluc- 
tuations on scales shorter (longer) than b/A, where A represents an ultraviolet cut-off and 
< b < 1. Applying the parameterization 
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Q = UQ > U-\ Q > =T > afT-\ U = RT < , 
where the rotation R defines the saddle-point (|55p, we obtain 



S = J str [D ((<9Q>) 2 + 4Q><9Q> + Q>] 2 ) + AieQ^afU 



(El) 



(E2) 



where $ = U dU. Setting T> = e , and expanding Q> to quadratic order in in the 
generators of rotations W, the action separates into three contributions: 

S = Ss + Ssf + <Sf, 



where, defining Q< = Ua^ h U 1 and Q < = U V ph 



3 ! 



5s 



TTUr, 



str 



£> (0Q<)* + 4zeaf Q< 



(E3) 



(E4) 



5 SF = — 7ru n J str dUw, dW] ■ $ + + ) 2 W^ 2 - 9iy ■ $ + iy($a : 

+ ie (WafQ < + W 2 afQ 



Ph\2 
3 J 



s F 



Dstr 



{dwf + %w 2 



(E5) 
(E6) 



We now integrate over the fast fluctuations, at the one-loop level, to obtain a new effective 
action, of the form S = Ss + (5sf)f- Taking the energy cut-off A/b to be far above E c , it is 
sufficient to neglect the linear terms in W and to integrate over the whole range of energies 
[A/6, oo) with a constant R. We find 



(Ssf)f 



7TZA, 



dr DU(r,r)str(dQ < f 



where n(r, r') represents the diffusion propagator 

2-KU n D (-d 2 + fpj Il(r, r') = 6(r - r'). 
Altogether, applying the rescaling, at one-loop we obtain the renormalized action 



S' 



str 



D cS (dQ) 2 + AieafQ 



(E7) 



(E8) 



(E9) 



where D c g = D[l — 211(0, 0)] denotes the renormalized diffusion constant. From this result, 
we see that the bare diffusion constant is subject to a standard weak localization correc- 
tion WM, albeit derived from purely within the particle/hole sector. 

As a result of the renormalization procedure, no new terms are generated in the effec- 
tive action. Instead, we obtain the usual kinetic term but with a renormalized diffusion 
coefficient. However, as a consequence of the renormalization, q • o"3 h = R~ 1 a^ h R no longer 
represents the saddle-point of the theory. Accordingly, it is necessary to recalculate the 
saddle-point solution in the presence of the renormalized diffusion constant. The result is a 
corresponding renormalization of the minigap edge discussed in the text. 
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